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1.0 Introduction 


CARE III (Computer-Aided Reliability Estimation, version 
three) is a computer program designed to help estimate the re- 
liability of complex, redundant systems. Although the program 
can model a wide variety of redundant structures, it was developed 
specifically for fault-tolerant avionics systems - systems dis- 
tinguished by the need for extremely reliable performance since a 
system failure could well result in the loss of human life. 

It IS usually relatively easy to design enough redundancy 
into a system to reduce to acceptably small levels the probability 
that It fails due to inadequate resources. The dominant cause of 
failure in ultra-reliable systems thus tends to be due not to the 
exhaustion of resources but rather to the failure to detect and 
isolate a malfunctioning element before it has caused the system 
to take an erroneous action. Such failures are called coverage 
failures. CARE III differs from its predecessors in, among other 
things, the attention given to coverage failure mechanisms. 

The first CARE program, developed at the Jet Propulsion 
Laboratory in 1971, provided an aid for estimating the reliability 
of systems consisting of a combination of any of several standard 
configurations (e.g. standby-replacement configurations, triple- 
modular redundant configurations, etc.). CARE II was subsequently 
developed by Raytheon, under contract to the NASA Langley Research 
Center, in 1974. It substantially generalized the class of redun- 
dant configurations that could be accommodated, and included a 
coverage model to determine the various coverage probabilities as 
a function of the applicable fault recovery mechanisms (detection 
delay, diagnostic scheduling interval, isolation and recovery 
delay , etc . ) . 
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CARE III further generalizes the class of system structures that 
can be modeled and greatly expands the coverage model to take into account 
such effects as intermittent and transient faults, latent faults, error 
propagation, etc. In order to accomplish this, it was necessary to 
depart substantially from the approaches taken in previous reliability 
modeling efforts. The mathematical reliability and coverage models 
implemented in CARE III are described in the CARE III PHASE II Report, 
Mathematical Description.* The following section, taken from that 
report, is included here in order to put into context the ensuing pro- 
gram maintenance discussion. 


*J. J. Stiffler and L. A. Bryant, CARE III PHASE II Report, 
Mathematical Descic.mtioQ. NASA CR-3566, 1982. 
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2 . 0 Mathematical Model Description 


The following paragraphs describe in detail the 
mathematical model as it is implemented in CARE III . The 
system to be modeled is assumed to consist of some number (up 
to 70) stages with each stage composed of one or more identical 
interchangeable elements or modules. The modules in each stage 
are subject to up to five user-defined categories of faults. A 
fault IS characterized in terms of its rate of occurrence and in 
terms of its coverage model parameters. Fault occurrence rates 
are constrained to be of the form toAt^ ^ (i.e., fault distribu- 
tions are constrained to be Weibull) with w and X user defined. 

The user can also specify up to five sets of coverage model para- 
meters (a, 6 , p(t), e(t), 6 (t) , C, P. , P„) ; each such set defines 

A D 

a fault type. (Thus, for example, it is possible to define a 
permanent fault type, a = 0 ; a transient type, a ^ 0, P = 0 ; and 
an intermittent a 7 ^ 0 , B ^ 0 ; each having its own characteristics 
with regard to detectability, error-propagation, etc.) Fault 
category x^ then refers to a fault that can affect any module in 
stage x; it is characterized by the parameters X , u) , and j 
with j a fault- type designator. 

In addition, the user must specify the number of modules n 
initially available at each stage, the minimum number m^ needed 
for that stage to function properly, the various combinations of 
stage failures that constitute a system failure, and the proba- 
bilities b (V ,v ) that a specific module in stage x forms a 
critical pair with a specific module in stage y given that v 
stage-x modules and stage-y modules are known to have failed 
and are therefore no longer being used.* 


♦These last two tasks are both accomplished with relative ease 
through a CARE III user interface incorporating a program called 
FTREE developed by Boeing Aircraft Co. and described in the CARE 
III User’s Manual. 
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On the basis of this user-supplied inforination, CARE III 
then determines the system unreliability using the equation 

R(t) = 1 - R(t) = 2 Q^(t) + 2 

JleL £eL 


with L the set of module failure combinations that would leave 
the system operational in the absence of a coverage failure, L the 
complementary set, P^*(t) the probability that the system would 
be in state Z at time t in the absence of a coverage failure, and 


ft 

f 

i If % 

- 

II 

( 

o 

T 

C (x[£-r )P *(x)(n -£ +1) A (x) 

JO 1 

I V . 

y . — y £-e y y y . 

- y D 
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+ A' (t|^)P^* (t) + a’ (t|£)Pj^*(t) 


dr 


wi th 


Uj^(t) = A' (t-l£) + a' (t|£) 


( 2 ) 

A .„(t) = / C (T|«-e ) (n -9. +1)A (x) 

]_£ il. y. y y y y 

y^ey ^ 


and 2 = —~^y unit vector denoting a stage-y module. 

( 2 ) 

Vi^(t) and Xj^(t) are defined in terms of functions of the form 


S, 


?2 (t)Pi ( t-T) dx 


with Pj, (t) a measure of the rate at which a certain class of faults 
occurs and P 2 (t) a function of the interval x between that occur- 
rence and the entry of the fault into a particular coverage-model 
state. Since, typically, faults occur at rates no greater than 
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one fault every several thousand hours, and since coverage-state 
time constants are usually of the order of fractions of seconds 
and rarely exceed a few minutes in duration, Pj^(t) is a much more 
slowly varying function of time than is • Thus, to a very 

good approximation 

2 

Pj^(t-i) sj a(t) + Tb(t) + T c(t) 

over the range of t for which Pj^(t) is not negligibly small, with 
a(t), b(t), and c(t) suitably defined- This approximation is used 
in CARE III with a(t), b(t), and c(t) defined to make the approxi- 
mation exact at the two end points and at the midpoint of the range 
of interest of t. The major advantage of introducing this approxi- 
mation is that, with it, 

^ (i) P^ (t-T) dx Si a(t)m2^(t) + b(t)m2^(t) + c(t)m2^(t) 

with 

m2^(t) T^P2(x)dx 

Thus, the convolution can be separated into two parts, one part 
depending only on the reliability-model function Pj^(t) and the 
other involving only the first three moments of the coverage-model 
function P 2 (t). Moreover, these moments need be evaluated only at 
those points of time t relevant to the reliability model. This 
significantly simplifies the interface between the coverage and 
reliability models. 

With these preliminaries, the reliability model functions 
used in CARE III are itemized in Table 3 and the coverage model 
functions in Table 2. 
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LIST OF SYMBOLS 


Symbol 

Module (s) 

Definition 

A(t|M 

CARE 3 

See Table 1 . 

A' (t|£) 

CARE3 

See Table 1. 

a (t) 
X 

CARE 3 

See Table 1 . 

a (t) 

CARE 3 

See Table 1 . 

a'(t|£) 

CARE 3 

See Table 1 . 

B (t|A) 

x.,y. 

CARE 3 

See Table 1 . 

b (V ,v ) 

x,y x' y 

• 

CARE IN 
CARE 3 

Probability that a stage-x fault and a 
stage-y fault form a critical pair given 
V known stage-x faults and v known 
stage-y faults . ^ 

c 

CARE IN 
COVRGE 

Probability that a propagated error is 
detected before it causes the system to 
malfunction . 

C (t|A) 

CARE 3 

See Table 1. 


CARE 3 

See Table 1 . 

•’x. 

1 

COVRGE 

Probability that a category-x^ error has 
not been detected t time units after its 
generation. 

e (t) 

COVRGE 

Probability that a category-x. error has 
not yet propagated t time units after 
its generation. 

Hjj(tlx^) 

X = L,B,B 

CARE 3 

See Table 1 . 

hx<t|Xi) 

X = DPT,F 

CARE 3 

See Table 1 . 

hnp(t|x^,yj) 

CARE 3 

See Table 1 . 
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Symbol 

Module (s) 

Definition 

i 

CARE3 

A vector, the components of which indi- 
cate the number of faulty elements in 

each stage (^ =...£, Jl . 

X y 


CARE 3 

( • • • • AjjP « • • • ) 

''x 

CARE 3 

Number of stage-x modules that have 
experienced some fault (X, = Z 1 ) 

X j X • 

i X 

s 

CARE 3 

Number of stage-x modules that have 
experienced category-x^ faults. 

m 

X 

CARE IN 
CARE 3 

Minimum nvimber of units needed for 
stage x to be operational . 

“x 

• 

CARE IN 
CARE 3 

Number of functional stage-x units 
available at time t = 0 . 


CARE 3 

State of operational system after ex- 
actly 1 failures . 

P(t) 

CARE 3 

Element survival probability. 

P'(t) 

CARE 3 

Rate of change of P(t) . 

p,(t) 

CARE 3 

Probability that a system has sustained 
exactly ^ failures and is still opera- 
tional at time t. 

p,‘(t) 

CARE 3 

Probability that a system has sustained 
exactly £ failures by time t. 

p(M^,tly 

CARE 3 

See Table 1. 

'■a 

CARE IN 
COVRGE 

Probability that a fault detected in the 
active state is diagnosed as permanent. 

*■8 

CARE IN 
COVRGE 

Probability that a fault detected in the 
benign state is diagnosed as permanent . 


CARE 3 

State of system that has failed after 
exactly ^ element failures . 
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Symbol 

Module (s) 

Definition 

Ql(t) 

CARE 3 

Probability that a system has sustained 
exactly ^ failures and then malfunctioned 
by time t. 

r (t) 

COVRGE 

Probability that a category-x. fault has 
not resulted in an error t time units 
after it has become active. 

R(t) 

CARE3 

Reliability; probability that the system 
is operational at time t. 

■^Ij <"> 


Transition rate between states S. and S.. 

1 D 

R^lt) 

CARE 3 

See Table 1. 

X 

• 

CARE3 

See Table 1. 

S. 

X 


State i in a Markov chain. 

Sj^(t) 


State i occupancy probability at time t. 

S^' (t) 


Rate of change of S^(t) . 

X 

CARE IN 
CARE 3 

Stage index. 

a 

X. 

X 

CARE IN 
COVRGE 

Rate of transition from active to benign 
fault states . 


CARE IN 
COVRGE 

Rate of transition from benign to active 
fault states. 


6 (t)/d (t) COVRGE 

Jv • A • 

X X 


Rate of detection of type-x. faults t 
time units after they become active. 


e (t)/e (t) COVRGE 

X f X • 

X X 


Rate at which an error caused by a type- 
fault is propagated following its gen- 
eration at time t = 0 . 
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Symbol 

Module (s) 

Definition 


CARE 3 

J X^(n)dn 



Rate of occurrence of failures that take 
the system from state P. to state P. . 

1 L 



Rate of occurrence of failures that take 
the system from state P^ to Q 




X,(t) 


Transition rate out of state P^. 

X,*(t) 


2-j x,.*(t) 

'« 

CARE IN 

Rate of occurrence of category -x. faults 

X 

COVRGE 
CARE 3 

u -1 

(Xx^^>=“x ^x ^ )• 

i i i 

y 

CARE3 

A vector, the components of which indi- 
cate the current number of latent faults 

in each stage (p=. . .y,vi, . . . 

X y 

P„ (t)/r (t) 

Xi x^ 

COVRGE 

Error generation rate of a type-x^ fault 
t time units after it becomes active. 

Vj(t) 

CARE 3 

Transition rate from system state P^j^ to 
system state Q, ^ ^ ^ _ 

(1) 

“4<« 

(2) 


Transition rate from state Pj^ to state 
due to single-fault coverage failures. 



Transition rate from state P^j^ to state 
due to double-fault coverage failures, 
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Reliability Model Functions 


FUNCTION 

H^(t) 
a (t) 

a^(t) 


MATHEMAT I CAL EXPRESSION DEFINITION 


I -A„.(t) 

1 e 


PERMANENT 


TRANSIENT 

n R (t) 
1 1 



PERMANENT 

l-R^(t) 


HL(t |x^) 

TRANSIENT 

^^a ct) 

> 

1 ^i 



(permanent) 


PROBABILITY THAT A GIVEN STAGE 
X MODULE HAS NOT EXPERIENCED A 
CATEGORY-Xj^ FAULT BY TIME t 

RELIABILITY OF A STAGE x MODULE 


PROBABILITY THAT A GIVEN STAGE 
X MODULE HAS A CATEGORY"X^ 
LATENT PERMANENT (TRANSIENT) 
FAULT AT TIME t GIVEN THAT IT 
HAS (not) experienced A PER- 
MANENT OR LEAKY TRANSIENT FAULT 
BY TIME t 

f 

PROBABILITY THAT A GIVEN STAGE 
X MODULE HAS A LATENT (PERMAN- 
ENT) FAULT AT TIME t GIVEN THAT 
IT HAS EXPERIENCED SOME PER- 
MANENT FAULT BY TIME t 




DEFINITION 


PROBABILITY THAT A SUBSYSTEM 
CONTAINS STAGE x LATENT 
PERMANENT FAULTS GIVEN THAT 
IT HAS 'jj STAGE x PERMANENT 
FAULTS 

EXPECTED NUMBER OF 
CRITICAL FAULTS AT TIME t 
GIVEN £ PERMANENT FAULTS 



Table 1 (Continued) 


Reliability Model Functions 


FUNCTION MATHEMATICAL EXPRESSION 


DEFINITION 


C(x^,y^) 
(CONT. ) 



(t|i) 
^ j 



i )u 
x' 


(t) 


= TRANSIENT 
Yj = PERMANENT 


^i'^j " TRANSIENT 

X y 


(n ) (n -S. -1) 

X X X X 



y 


X 


y 


y 



= TRANSIENT 

X = y 


,a -V ) P (y , t U ) 

' y y ' x' ' X 


P(Py,t|«.y) 

a (t) 

X. 

^x a W ""i = permanent 

X 

= transient 


^i Hj^(t|x^) °x^,y 


EXPECTED NUMBER OF x^y- 
CRITICAL FAULTS. GIVEN 1 

permanent faults, that would 

BE CREATED AS THE RESULT OF 
A STAGE y FAULT AT TIME t 


PROBABILITY THAT A CATEGORY y^ 
FAULT WOULD PRODUCE A SYSTEM 
FAILURE AT TIME t GIVEN L 
FAULTS AT TIME f 



Table 1 (Continued) 


Reliability Model Functions 


FUNCTION MATHEMATICAL EXPRESSION 


DEFINITION 


a' (t|_a) 


a (t[^) 


H^(tlx.) 


H3(tlXi) 


B (t|e,) 


S ^j^hp(t|Xi) 

^i 1-R (t) 

X 

PERMANENT 

^i 

TRANSIENT 


(t 1 x^)M° (t I +bj^ (t 1 x^)M^(t 1 x^) 
+ Cj^(t|x^)M^(t|x^) 




RATE WHICH SYSTEMS HAVING i 
FAULTS FAIL AT TIME t DUE TO 
CRITICAL FAULT CONDITIONS 

RATE AT WHICH SYSTEMS HAVING 
1 FAULTS FAIL AT TIME t DUE TO 
ERROR PROPAGATION 


PROBABILITY OF A LATENT CATEGORY 
FAULT AT TIME t 


PROBABILITY OF A NON-BEN IGN 
LATENT FAULT AT TIME t 



Table 1 (Continued) 


Reliability Model Functions 


FUNCTION MATHEMATICAL EXPRESSION DEFINITION 


HgCtlx.) 


ag(t| x^)M° (t|x^)+bg(t|x^)Mg(t |x^) 



PROBABILITY OF A BENIGN 
LATENT FAULT AT TIME t 




DF 


(tlx^, 




(t|x^) 


aDp(t|Xi)m°p(t|Xi)+b^p(t|Xi)mJp(t|Xi) 


+ Cj^p(t|x.)mj^p(t|x^) 


PROBABILITY THAT A CATEGORY 
TRANSIENT FAULT IS DETECTED AS 
PERMANENT 


+ bj3p(t|x.,y^)m^p(t|x.,y^). 

+ Cj^p(t|Xi,yj)m2F(t|Xi,y 1 


ap (t I x^)mp (t I x^)+bp (t |x^)mp (t ] x^) 


+ Cp(t|x^)mp(tlx^) 


RATE AT WHICH AN x^y^-CRITICAL 
FAULT CAUSES SYSTEM FAILURE 


RATE OF ERROR PROPAGATION 
FAILURE DUE TO A CATEGORY 
FAULT 



FUNCTION 


Reliability Model Functiong. 

MATHEMATICAL EXPRESSIOM DEFINITION 


ajj(t |Xj^) 
aj;(t|Xi,y.) 



all X 


f (t) 




X = 


X. = 
1 


X (t) R (t) 

X. X 

1 

X. = 
1 


HB(tlXi)Xy 


y. = 

j 


HB(t|x.) X 

j 

y. = 

3 


DP,L,B,B,F 

TRANSIENT 

X = L,B,B,F 
NON-TRANSIENT 

(t) X = DF 

TRANSIENT 

(t)r^(t) X = DF 

NON-TRANSIENT 


VJEIGHT FUNCTIONS USED IN 
CONVOLUTIONAL APPROXIMATION 
(CF. EQUATION 25) 


f(t) - f(0) 


t = £.t 


(3k^-l)f(t) - 


+ (k^+l)f(0) 


(k -l)t 


t = kAtj., < tj^, k odd 


3f(t) - 4f(t/2) + f(0) 


t = kAt , < t , k even 
r n 


3f(t) - 4f(t-t /2) + f(t-t ) 
n n 


t 

n 


t > t 


n 


V 



Table 1 (Contxnued) 
Reliability Model Functions 


FUNCTION MATHEMATICAL EXPRESSION 


DEFINITION 


f(t) - 


t = At 






+ f(0) 


[t^ - (4t^)^]/2 


o^(t|x.) 


Cj((t|x^,y.) 


t = kAt < t 
r n 

k odd 


f(t) - 2f(t/2) + f(0) 
t ^/2 


t = kAt < t 
r n 

k even 


f(t) - 2f(t-t /2) + f(t-t ) 
n n_ 

t ^/2 
n 


t > t 


n 


t 

n 


min t = n t , n even, such that P„(t) < 9 

IT X 

or P (t) 9 with 9 a user-defined threshold 



Table 2a 


Single-Fault Model Equations 


FUNCTION 


MATHEMATICAL EXPRESSION^ 


DEFINITION 


(t) 


ae 


-et 


f (x)d(T)dT 

Jo 


TIMES THE PROBABILITY 
INTENSITY OF RE-ENTERING 
STATE A EXACTLY t TIME 
UNITS AFTER THE PREVIOUS 
ENTRY 


P^(t) 


e”“^r(t)d(t) + B \ 4>(t-T)P (x)di 

a 


s: 


PROBABILITY OF BEING IN 

STATE A AT TIME t WHEN 

P = P =1 
A B 


Pb(t) 


j 0 


(t) + F \ ^ (t-T)Pj^(T)dT 


PROBABILITY OF BEING IN 
STATE B AT TIME WHEN 

P = P = 1 
A B 




s: 


-at 


e P (T)d(T)e(t-T)dT + 8\ -t(t-x)P (x)dx 




PROBABILITY OF BEING IN 
STATE Ac OR B. AT TIME t 
WHEN " ^B " ^ 


t HERE IS A MEASURE OF THE TIME SINCE THE ENTRY INTO STATE A. 



TaDle 2a (Continued) 


Single-Fault Model Equations 


FUNCTION MATHEMATICAL EXPRESSION* 


DEFINITION 


Pe<t) 


Pg(t) 


00 


Pf (t) 


4^(t) 


e"“^p(t)d(t) + p 


e (5( t) r (t) + B 



4>(t--i)p^(T)dT 


4i (t-T) p-(t) di 


INTENSITY OF ENTRY INTO 
STATE Ag AT TIME t WHEN 

P = P =1 
A B 

INTENSITY OF ENTRY INTO 
STATE Ap FROM STATE A AT 
TIME t WHEN ^ 


(1-C) 



p (T)e(t-T)dT 


INTENSITY OF ENTRY INTO 
STATE F AT TIME t WHEN 


^ P^(t) c (t-T)^— 


-(a+P) (t-T) 


a+ P 


) 


dx + p-(t) 
e 


INTENSITY OF ENTRY INTO 
STATE Ajj AT TIME t FOR 
THE FIRST TIME 




t HERE IS A MEASURE OF THE TIME SINCE THE ENTRY INTO STATE A. 



Table 2a (Continued) 


Single-Fault Model Equations 


FUNCTION MATHEMATICAL EXPRESSION* 


DEFINITION 




X B 



Pg ( t) ( 1 


- (a+B) (t-T) . 

e ) c 


( t-T ) dx 


INTENSITY OF ENTRY INTO 
STATE Bq at time t FOR 
THE FIRST TIME 



Yg(x)e"^^^ ^^dx 


PROBABILITY OF HAVING 
ENTERED STATE B^ FOR THE 
FIRST TIME AND THEN RE- 
MAINING IN THE BENIGN 
STATE UNTIL TIME t 


P 


dp 


(t) 


A 



(x)dx + 



Yjj(T)dl 


PROBABILITY THAT A FAULT 
HAS BEEN DIAGNOSED AS 
PERMANENT BY TIME t 


F^(t) 


« 


r 

Jo 


F„(t) +\ [ (l-P^)H'^ (t-x) + (1-Pg) BxB^t-x) ]F^(x)dx 


A' A 


X 


FUNCTION RELATING PROB- 
ABILITIES AND INTENSITIES 
DERIVED WHEN P^ = Pg = 1 TO 
THOSE SAME QUANTITIES WHEN 
^A ^ are ARBITRARY 


t HERE IS A MEASURE OF THE TIME SINCE THE ENTRY INTO STATE A. 



Table 2a (Continued) 



FUNCTION MATHEMATICAL EXPRESSION* DEFINITION 


Pg(t) F^it) with F^(t) = Pj^(t) + Xg(t) PROBABILITY OF BEING IN 

STATE B AT TIME t 

PB(t) F^(t) with F^(t) = P^(t) + P^(t) PROBABILITY OF BEING IN 

A NON-BEN IGN STATE AT 
TIME t 

Pb(t) + Xg(t) 

+ P (t) + P (t) 

a e PROBABILITY OF A LATENT 

PERMANENT FAULTS FAULT OR UNDETECTED ERROR 

AT TIME t 

P^(t) + P .(t) 
a > e 

TRANSIENT FAULTS 

^DP^^^ ^X^^^ PROBABILITY THAT A FAULT 

HAS BEEN DIAGNOSED AS 
PERMANENT BY TIME t 

* t HERE IS A MEASURE OF THE TIME SINCE THE ENTRY INTO STATE A. 




Table 2b 



Equations 


DEFINITION 

TRANSITION RATE FROM 
STATE AjB3^ TO STATE F 


TRANSITION RATE FROM 
STATE A.B. TO STATE 

J 


INTENSITY OF ENTRY INTO 
STATE F t TIMEUNITS AFTER 
ENTRY INTO STATE B^B2 


INTENSITY OF RE-ENTRY 
INTO STATE t TIME 

UNITS AFTER A PREVIOUS 
ENTRY 



Table 2b (Continued) 


Double-Fault Model Equations 


FUNCTION MATHEMATICAL EXPRESSION 'DEFINITION 


P3(t) f^it) + \ C3(t-T)P3(T)dT INTENSITY OF ENTRY INTO 

STATE B 1 B 2 t TIME UNITS 
AFTER ENTRY INTO STATE 

A2B1 

Pop(t) c^^It) + \ c^(t-x)p3(T)dT INTENSITY OF ENTRY INTO 

STATE F t TIME UNITS 
AFTER ENTRY INTO STATE 



3.0 CARE III Program Modules 


CARE III consists of three interdependent software modules written in 
FORTRAN EXTENDED Version 4: 

1) CAREIN - input processor, 

2) COVRGE - coverage model , and 

3) CARE3 - reliability model (see Figures 3.0.1 and 3.0.2). 

The following sections contain a description sheet for each of the main 
programs, subroutines, and functions which comprise each of the modules listed 
above. Each routine description sheet consists of the files, parameters, 
arrays, common areas, externals, and inline functions required by the 
routine. Where appropriate, the mathematical functions are included with the 
routine that solves them. Afunctional flowchart is included for each module, 
and a description of the files, enumerated in the program card of each main 
program, is also included. 

The following are the execution field lengths for each program com- 
prising CARE III. 


MODULE 

PROGRAM 

FIELD LENGTH (Octal) 

CAREIN 

CAREIN 

154000 

COVRGE 

COVRGE 

163700 


CVGPLT 

127600 

CARE3 

CARE3 

157200 


RELPLT 

076000 


23 



CARE III STRUCTURE 


SUCCESS CRITERIA 



USER ' 
INPUTS^ 


SYSTEM 


STRUCTURE; 
FAULT PARAMETERS 


STATE 
TRANSITION 
RATES 


SYSTEM 
RELIABILITY 


Figure 3»0.1 





























3.1 CARE IN Module 


The CAREIN software module consists of the main program CAREIN 
and several subroutines and functions. A major portion of this 
module consists of subroutine FTREE and its related routines de- 
veloped by Boeing Aircraft Company and described in the CARE III 
PHASE II Report, User's Manual.* 

The CAREIN module processes user input contained in file 
CREIN and then generates files for modules COVRGE and CARE3. (See 
the CARE III User's Manual for a description of the user input to 
CAREIN.) 


*L. A. Bryant and J. J. Stiffler, CARE III PHASE II Report, 
User's Manual, NASA CR-165864, September 1982. 
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3.1.1 CARE IN Main Program 

PROGF<AM CAREIN( INPUT T OUTPUT fTAPE5--=INPUT»TAPE6 = 0UTPUT>CREINj 

1 TAPE7=CREIN»CREOUT>TAPE8=CREOUT»TAPE10=82»FT11F»TAPE11=FT11F» 

2 BXYFL = 0»TAPE14 = BXYFL» FT15F»TAPE15 = FT15F,FT16F» TAPE16=FT16Fy 

3 FT25Ff TAPE25=FT25FfFT26Ff TAPE26=FT26F»C0VIN=0rTAPE3=CQVIN» 

4 RELIN=0»TAPE4=RELIN) 

C 

m ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ 4 ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ 

^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ *. ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ <T* ^ * * * * ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ 


c* * 
C* title: computer-aided reliability ESTIMATION? VERSION THREE * 
C* acronym: CARE III -1^ 
C* MODEL design: DR. J.J.STIFFLER ♦ 
C* SOFTWARE DESIGN AND IMPLEMENTATION: LYNN A. BRYANT * 
C* date: MARCH 23?1981 ♦ 
C* DATE edited: 4-22-81 * 
C* EDITED by: LYNN A. BRYANT * 
Z% ♦ 

c* purpose: care hi is a general-purpose reliability * 
c* estimation tool for fault-tolerant avionics * 
c* ' systems. * 
c* * 


^ "4f ^ ^ ^ ^ ^ ^ ^ 4 ^ 4 ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ 4 ^ 4 ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ 4 ^ ^ ^ ^ ^ ^ ^ 4 ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ 4 ^ 

^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p o* ^p 

c 

C FT15F IS the system MINTERM FILE. 

C FT25F IS THE CRITICAL PAIR MINTERM FILE. 

C 

COMMON// P(2000) ?T(50) ? SUM (2000) ? PR ( 70 ) ? LS ? LTT ? S AV ? MSENS ? I D ? 

1 QP(70) ?PQ( 70) ? BOG? BIG? DELTA ? IC( 2000) ? IN(2000) ? IP( 10000) ? IR(2000) ? 

2 IDEP(700) ?IDUM(4300) ?CINDBG 
LOGICAL CINDBG?FATAL?EOFFLG 

c* note: 'idum' is a place-holder for use in subroutine ctrlprs when 

C* 'IC' IS EQUIVALENCED to 'KNTMTS" dimensioned to (10?10?2)0). 
c 


EXTERNALS 

TYPE 

ARGS 

bufblk 

' 

6 

CRTLPRS 


4 

EOF 

REAL 

1 

FNCK 

REAL 

2 

FTREE 


1 

SPLIT 


6 

SUBRUN 


6 


namelists 

FLTCAT 

flttyp 

RNTIME 

STAGES 
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3.1.2 CAREIN Files 



File Name 

Unit Name 

Description 

i! 

COVIN=i? , 

TAPE 3 = COVIN, 

coverage definition input file 
generated by this program for 
COVRGE program. 

2! 

RELIN=j2f, 

TAPE4 = RELIN, 

system definition input file 
generated by this program for 
CARES program. 

3. 

INPUT, 

TAPES = INPUT 


4. 

OUTPUT, 

TAPES = OUTPUT 


5. 

CREIN, 

TAPE? = CREIN 


6. 

CREOUT, 

TAPES = CREOUT 


7. 


TAPEIO, 

TAPEll converted for READ* . 

8. 

FTllF, 

TAPEll = FTllF, 

FTREE input . 

9! 

BXYFL=j2f, 

TAPE 14 = BXYFL, 

binary critical pair results 
file passed to CARE3 program. 

lot 

FT15F, 

TAPE15 = FT15F, 

system minterm file passed to 
CARE 3 program. 

11. 

FT16F, 

TAPE16 = FT16F, 

fault tree output from system 
fault tree definition. 

12. 

FT25F , 

TAPE25 = FT25F, 

critical pair minterm file. 

13. 

FT26F, 

TAPE26 = FT26F, 

fault tree output from critical 


pair fault tree definition (s) . 


^Must be saved and passed to other programs . 
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3.1.4 CARE IN Subroutines 


^i/ ^ «Jf 1^ ^U ^U ^U ^b ^ ^U ^b 

^n ^p ^p ^p ^p ^p ^p ^n ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p 

C 

C* SUBROUTINE BINOM 

c 

C* THIS SECTION CALCULATES THE PROBABILITY OF 

C* M AT LEAST N. H HAS A MAXIMUM VALUE OF 20. 

C* PROGRAM WILL HANDLE UP TO 10 N'S AND THEN 

C* DOES A FLIPFLOP TO CALCULATE AN N > 10. 

C 

c 

C* VARIABLE DEFINITIONS 

C 


c* 

Q1 

20 

PROBABILITY 


c* 

Q2 

20 

COMPLEMENT 

OF PROBABILITY 

c* 

ST 

1 

CALCULATED 

PROBABILITY FOR EACH MINTERM 

c* 

C 

20 

COUNTER AND 

CONTROLLER 

c* 

NT 

1 

NUMBER OF TERMS(M) 

u 

c* 

JJ 

1 

MINIMUM NUMBER OF SUCCESS(N) 


C 

/*% ^ ^b ^b ^b ^b ^b ^b ^b ^b ^b ^b ^b ^b ^b ^b ^b ^b ^ia ^b ^b ^b ^b ^b ^b ^b ^b ^b ^b ^b ^b ^b ^b ^b ^b ^b ^b 

^P ^P ^P ^P ^P ^P ^P ^P ^P ^P ^P ^P ^P ^P ^P ^P ^P ^P ^P ^P ^P ^P ^P ^P ^P ^P ^P ^P ^P ^P ^P ^P ^P ^P ^P ^P ^P ^P ^P ^P ^P ^P ^P ^P ^P ^p ^P 

SUBROUTINE BINOM <NT»JJ»QfST> 

REAL Q(20) iQl <20) »Q2(20) »ST 


EXTERNALS 

TYPE 

ARGS 

ec 

CALC 


D 
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c 

C* SUBROUTINE CALC 

C* CALLED BY BINOM 

C 

C* DOES ACTUAL CALCULATION FOR *AT LEAST* CONDITION* 

C 

j|c *1^ i|c ^|c ^|c }^ ^C ^|c ^|c 

SUBROUTINE CALC ( Q1 » Q2 » ST f C » NT ) 

REAL 01 (20) jQ2(20) rSSfSfST 
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SUBROUTINE CRTLPRS ( N i M » NOP » LC ) 

C 

COMMON// P<2000) fT(50) »SUM(2000) »PR(70) »LS»LTT»SAV»MSENS»ID» 

1 QP(70) >PQ(70) »BOG»BIG»DELTAf IC<2000) , IN(2000) , IP( 10000) » IRC2000) r 

2 IDEP(700) »IDUM(4300) rCINDBG 
LOGICAL CINHBG»FATAL»EOFFLG 

DIMENSION KNTMTS( 10» 10»210) »MNTRMV(70> ,MTSTGS(70) 

EQUIVALENCE < IC < 1 ) » KNTMTS ( 1 f 1 f 1 ) ) 

C 

COMMON/BXYCOM/ BXYAR < 10 r 10 » 55 ) » IJSTGIN ( 210 ) » LCNDVEC ( 20 ) » IBREC » 

1 KFSTGfKPLSTG(70> 

DIMENSION BC0MAR(5802) 

EQUIVALENCE < BX YAR ( 1 » 1 » 1 ) » BCOMAR < 1 ) ) 

DIMENSION N(70) »M(70) »N0P(5»70> fLC(70) 


EXTERNALS 

BUFBLK 

TYPE 

ARGS 

6 


EOF 

REAL 

1 


FNCK 

REAL 

2 

- 

INLINE FUNCTIONS 

TYPE 

ARGS 


FLOAT 

REAL 

1 

INTRIN 

MINO 

INTEGER 

0 

INTRIN 
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SUBROUTINE DFIELD ( ARRAY » 1 D I H » I TAPE r LN ) 
C* SUBROUTINE DFIELD READS FREE FIELD FORMAT 
C 

DIMENSION ARRAY(50) >GATE(7) 


EXTERNALS 

EOF 


TYPE ARCS 
REAL 1 
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SUBROUTINE FTREE(IUNIT) 

^|c > 1 ^ ^ 

C* title: simplified computer evaluation of FAULT TREES 

c* acronym: sceff 

c* author: c. r. stander 

C* date: OCT, 25f 1976 

C* DATE edited: 4-22-81 

C* EDITED by: LYNN BRYANT 

C 

C* purpose: SCEFT provides the analyst REQUIRED TO WORK WITH 

C* FAULT TREES A CAPABILITY TO OUICKLY AND SIMPLY 

C* EVALUATE THE PROBABILITIES OF ALL EVENTS ON THE 

C* TREE, 


TYPES OF DEPENDENCIES THE PROGRAM UILL HANDLE ARE: 

INPUT EVENTS THAT OCCUR IN MULTIPLE PLACES IN THE TREE 
. COMMON TIME-SEQUENTIAL REDUNDANT PROBABILITIES 
(CONDITIONAL-AND) 

COMMON ACTIVE AND PASSIVE REDUNDANCY WITH OR WITHOUT 
SWITCHING DEVICESf AND 
EXCLUSIONS (EXCLUSIVE-OR) , 


EXTERNALS 

TYPE 

ARGS 

ALOG 

REAL 

1 LIBRARY 

DFIELD 


4 

EXP 

REAL 

1 LIBRARY 

MINTRMS 


4 

PROCIND 


0 

SORT 


2 

TTREE 


0 

INLINE FUNCTIONS 

TYPE 

ARGS 

ABS 

REAL 

1 INTRIN 
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n n n r> 


SUBROUTINE GNMNFL ( ZSUM , I UNI T » FRSTCM ) 

THIS SUBROUTINE WRITES VALID SYSTEM MINTERMS TO FILE ' FT15F ' ( IUNIT = 
15) OR VALID CRITICAL PAIR MINTERMS TO FILE ' FT25F ' ( IUNIT = 25 ) . 

DIMENSION IIPP(70) 

LOGICAL FRSTCMf CINDBG 

COMMON// P(2000) »T(50) f SUM (2000) » PR(70) » LS » LTT » SAV » MSENS » ID r 

1 QP(70) »PQ( 70) » BOG » BIG » DELTA »IC( 2000) » I N ( 2000 ) , I P ( 1 0000 ) » IR(2000) » 

2 IDEP(700) r IDUM(4300) f CINDBG 

C* note: 'IDUM' is a place-holder for USE IN SUBROUTINE CTRLPRS WHEN 
C* 'IC' IS EQUIVALENCED TO 'KNTMTS' DIMENSIONED TO (10rl0»210)» 

C 
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SUBROUTINE MI NTRMS ( NN » TSUM » I UN I T » FRSTCM ) 

C 

LOGICAL FRSTCM»CINDBG 

COMMON// F(2000) » T ( 50 ) r SUM ( 2000 ) i PR ( 70 ) j LS » LTT » SAM» hSENS f ID » 

1 QP(70) »PQ( 70) » BOG » BIG » DELTA » ICC 2000) » IN (2000) f IP (10000 )»IR( 2000) r 

2 inEP(700) rIDUM(4300) »CINDBG 

C* NOTEJ 'IDUM' IS A PLACE-HOLDER FOR USE IN SUBROUTINE CTRLPRS WHEN 
C* 'IC' IS EQUIVALENCED TO 'KNTMTS' DIMENSIONED TO (10»10>210)* 

C 


C» 

C* 

C* 

C* 

C 


THIS SECTION CALCULATES EACH MINTERM AS EACH 
MINTERM FAILS. MCOMB DETERMINES THE NUMBER OF MINTERMS 
THAT WILL BE EVALUATED BY THE SYSTEM. THE MAXIMUM NUMBER 
FOR MCOMB IS 13. 


EXTERNALS 

GNMNFL 


TYPE ARCS 
3 


TTREE 


0 
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SUBROUTINE PROCIND 

C* THIS ROUTINE PREPROCESSES INPUT DATA SO THAT SUBROUTINE DFIELD CAN 
C* READ THE INFORMATION USING FORTRAN READ*. NOTE: NON-NUMERIC DATA 
C* IS CONVERTED TO A NEGATIVE INTEGER REPRESENTING THE HOLLERITH 
C* CHARACTER» INSTEAD OF DELIMITING THE CHARACTER BY 'NOT EQUAL SIGNS'. 
C 


EXTERNALS 

EOF 


TYPE ARGS 
REAL 1 
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^U ft^ ^U ^« ^f ^« *^ ^U ^U ^> ^b ^b Wf •#« ^ ^b ^b ^* ^u ^ 

J^ J^ i^ ^h 

C 

C* SUBROUTINE SORT 

C* PUTS GATES IN ASCENDING ORDER 

C 

C* VARIABLE DEFINITION 

C 

C* N 1 NUMBER OF TERMS TO BE SORTED 

C* C 50 TERMS TO BE SORTED 

C 

^b ^b ^b ^b ^b ^b ^b ^b ^b ^b ata ^b ^U ^a ala ^b Ja ^b ^b ^b ^a ^a ^b ^a ala ^b ^b ^Ja ^b ^b ^b ^b ^b ^U ^a ^b .ta ^a ^a ^a ^a a^ 

^n ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p Ip ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p 

SUBROUTINE SORT (N»C) 

INTEGER C(50)iD 
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SUBROUTINE SPL I T < NKSTGF » NKSTGL t MXCPLS » MXSBRN . I SBRN » I SBRAR ) 

C 

C THIS SUBROUTINE SPLITS A RANGE OF NON-COUPLED STAGES INTO 'MXCPLS' OR 
C LESS SETS OF APPROXIMATELY EQUAL SIZE, THE FIRST AND LAST STAGES OF 
C EACH SET FOR EACH SUBRUN ARE THEN STORED IN ' ISBRAR ' ( I r ISBRN ) » WHERE 
C 'ISBRAR' <1 fISBRN) = FIRST STAGE IN THE SUBSET AND 
C 'ISBRAR' (2»ISBRN) = LAST STAGE IN THE SUBSET? 'ISBRN' IS THEN 
C INCREMENTED, 

C 

DIMENSION ISBRAR<2»MXSBRN) 
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SUBROUTINE SUBRUN < KSSETS »NSTGES » MXCPLS » hXSBRN » ISBRN » I SBRAR ) 

C 

C THIS SUBROUTINE DIVIDES 'NSTGES'f SOME OF WHICH ARE COUPLED» I.E, 

C THEY DEFINE CRITICAL-PAIRS OF FAULTS, INTO SUBRUNS. 

C 'KSSETS'(IrlSBSET) DEFINES ALL CRITICAL-PAIR FAULT SUBSETS. 

C ** EACH SUCH SUBSET MUST START AND END A NEW SUBRUN. ** 

C SUBROUTINE SPLIT IS USED TO DIVIDE THE RANGE OF NON-COUPLED STAGES 
C INTO SUBRUNS. 

C 

DIMENSION KSSETS(2,MXSBRN) , ISBRAR ( 2 , MXSBRN ) 


EXTERNALS 

SPLIT 


type arcs 
6 
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w ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ 'T^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ i^ 

C 

C* SUBROUTINE TTREE 

C* CALCULATES PROBABILITIES DEPENDING ON 

C* TYPE OF GATE(IR) DCBB78 

C 

(^ >fc ^Hc ^c )fc 

c 


c* 

c* 

Q 

20 

CONTAINS PROBABILITY FOR EACH INPUT GATE 
THAT FEEDS INTO *AT LEAST' GATE 

c* 

SUM 

2000 

HOLDS PROBABILITY OF EACH LOGIC GATE 

c* 

ST 

1 

RETURNS PROBABILITY FROM 'BINOM' 

c* 

SS 

1 

TEMP STORAGE FOR PROBABILITIES 


j|c 3jc 5(c jjc 3(C jjt 3(c 3j£ 3)£ 3)C 3jC )jC 3jC 3jC JjC 3jC 3jC 3jC 3jC 3jt 3jC ifC 3jC 3^C 3(C ijC JjC Jj{ 3jC 3(C 3jt JjC 3jC 3jC )fC 3({ 3jl 3jC IjC 3^ 3^ 3jC 3jC 3jC 3jC 3(% 3j[ 3|{ 3|C 3jC 3^C 3^ 3^ 3jC 3jC 3jC 3|C 3^ 3(( 3(( 3jt 

SUBROUTINE TTREE 
DIMENSION Q(20) 

COMMON// P ( 2000) »T(50) »SUM( 2000) »PR( 70) rLS»LTT»SAVfMSENS» IDr 

1 QP(70) tPQ<70) f bog » big » delta r IC( 2000) > in ( 2000) r IP (10000) » IR(2000) y 

2 IDEP(700) y IDUM(4300) yCINDBG 
LOGICAL CINDBG 

C» note: 'IDUM' is a place-holder for USE IN SUBROUTINE CTRLPRS WHEN 
C» 'IC' IS EQUIVALENCED to 'KNTMTS' DIMENSIONED TO (10yl0y210). 

C 


EXTERNALS 

BINOM 


TYPE ARGS 
4 
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3.2 COVRGE Module 


The COVRGE module consists of two main programs, COVRGE and CVGPLT, and 
several subroutines and functions. Program COVRGE computes the state transition 
rates required by program CARE3, (See the CARE III PHASE II Report, Mathematical 
Description* for a description of the mathematical model used in the coverage 
model implementation.) Program CVGPLT creates the plot file of the single- and 
double-fault functions generated in this module. 

For further information regarding the mathematical equations used in the 
coverage model see the Appendix entitled "Computational Complexity in Solving 
the 'Stiff Vol terra Equation." 


*J. J. Stiffler and L. A. Bryant, CARE III PHASE II Report, Mathematical 
Description . NASA CR-3566, 1982. 
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3.2.1 COVRGE and CVGPLT Main Programs 

PROGRAM COVRGE ( INPUT i OUTPUT » TAPE5=INPUT » T APE6»OOTPUT » DBUG » 

1 CVGMTS=Of TAPE7=CV6MTSf TAPE8=DBUG»SNGFL=0»TAPE9=SNGFL» 

2 DBLFL=0»TAPE10=DBLFL»C0VIN=0»TAPE3=C0VIN) 


^ ^ ^JU ^If ^ ^J/ 

I. TP ^p ^P ^p ^p ^p ^p ^p ^p ^p ^p ^p ^k ^P ^P ^k ^p J^k ^p ^p ^p ^k kp ^p ^p ^k 


c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 


title: computer-aided reliability estimation^ version three 

acronym: care hi 

MODEL design: DR. J.J.STIFFLER 

SOFTWARE DESIGN AND IMPLEMENTATION: LYNN A. BRYANT 
date: march 23fl981 
DATE edited: 5-02-81 
EDITED by: LYNN A. BRYANT 

purpose: care hi is a general-purpose reliability 

ESTIMATION TOOL FOR FAULT-TOLERANT AVIONICS 
SYSTEMS. 


♦ 

* 

* 

% 

% 

% 

% 

% 

% 

% 

% 

t 

* 


^ ^k ^k ^k ^k ^k ^k 

^p ^p ^P ^p ^p ^k ^k ^p 4 ^ ^k ^p ^p ^p Jp ^p ^p 4 ^ ^k 4 ^ ^P ^P ^P ^P 4 ^ ^P ^p ^P 


COMMON/CVGSTEP/ ITSTPS » MAXSTP» RELSTEP » TBASE > NTYPS » DBLDF » TRUNC » 
1 INMAXfNSTPMXrDIFCHNG 

DIMENSION RECK4) »REC2(3) 

EQUIVALENCE <ITSTPS»REC1 ( 1 ) ) , ( NTYPS » REC2 ( 1 ) > 


C 

COMMON/FLTYPCM/ ITYP » JTYP , CVPRNT r CVPLOT H AXSTYP 
DIMENSION REC3(3) 

EQUIVALENCE ( CVPRNT » REC3 < 1 ) ) 

C 

LOGICAL DELCDNfRHOCDN»EPSCDN 

COMMON/RATECOM/ ALPHA(5) » BETA ( 5 ) » DELT A ( 5 ) » RHO ( 5 ) » EPSLN ( 5 ) f 

1 DELCDN(5) »RH0CDN<5) »EPSCDN(5) r 

2 PACNS(5) »PBCNS(5) »CCNS(5) » 

3 TALPMX(5) »TBETMX(5) . TDELMX ( 5 ) , 1 RHOMX ( 5 ) » IEPSMX(5) 
DIMENSION REC4(55) 

EQUIVALENCE ( ALPHA < 1 )» REC4 ( 1 ) > 


COMMON/CVRGCOM/ CMSTO ( 5 » 65 » 5 ) » CMST 1 ( 5 » 65 > 5 ) » CMST2 ( 5 , 65 » 5 ) r 

1 CMDF0(5»5»65) f CMDF1(5»5»65) » CMDF2 ( 5 » 5 » 65 ) r 

2 TZER0ST(5»5) »T0DF(5»5) HRANS(5> 


?7 /P20 

C FOR OPTIMIZATION PURPOSES DURING THE BUFFER OUT OF /CVRGCOM/i- 
C DEFINE AND EQUIVALENCE 'CVRGAR' TO ENCOMPASS THE ENTIRE COMMON AREA. 
DIMENSION CVRGAR(9805) 

EQUIVALENCE < CMSTO ( 1 » 1 » 1 ) » CVRGAR ( 1 ) ) 

LOGICAL TRANS fCVPRNTfCVPLOTf FATAL »EOFFLG 


EXTERNALS 

TYPE 

ARGS 

ALOG 

REAL 

1 LIBRARY 

BUFBLK 


6 

DBLFLT 


0 

PRNTCVG 


0 

SNGFLT 


0 
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Figure 3. 2. 1.1 
CARE III SINGLE FAULT MODEL 

a = active to benign 
transition rate 

/3 = benign to active 

transition rate 

P(t) = error generation rate 

<(t) = error propagation rate 

5(t) = fault detection rate 

t = time from entry into 
active state 

T = time from entry into 
error state 



Figure 3.2.1 .2 
CARE III DOUBLE FAULT MODEL 
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PROGRAH CVGPLT < INPUT » OUTPUT »TAPE5= INPUT fTAPE6=qUTPUT»PLFILE» 

1 TAPE4=PLFILE»SNGFL=0»TAPE9=SNGFL»DBLFL=0,TAPE10=DBLFL,C0VIN=0» 

2 TAPE3=C0VIN) 

C 

*lf ^ ^ ^ 

^ Y T* ^P ^P ^p ^P ^P ^P ^P ^P ^P ^P ^P ^P ^P ^P ^P ^P ^P ^P ^P ^P ^P ^P ^P ^P ^P ^P ^P ^P ^P ^P ^P ^P ^P ^P ^P ^P ^P ^P ^P ^P ^P ^P ^P ^P ^P ^P ^p 


c* - * 
C* title: computer-aided reliability ESTIMATIONr VERSION THREE * 
c* acronym: care hi * 
C* MODEL design: DR. J.J.STIFFLER * 
C* SOFTWARE DESIGN AND IMPLEMENTATION: LYNN A. BRYANT * 


C* 

C* 

C* 

C* 

C* 

C* 

C* 

C* 


date: march 23»19S1 
DATE edited: 4-30-81 
EDITED by: LYNN A. BRYANT 

purpose: care hi is a general-purpose reliability 

ESTIMATION TOOL FOR FAULT-TOLERANT AVIONICS 
SYSTEMS. 


* 

* 

* 

* 

* 

* 

* 

* 


^|c 3^ ^|c ^|c ^|c 3|c }|c 

c 

C THIS PROGRAM PLOTS THE COVERAGE SINGLE AND DOUBLE-FAULT FUNCTIONS. 
C IT BUFFERS IN FROM FILE 'SNGFL' THE TWENTY-FIVE POSSIBLE SINGLE- 
C FAULT FUNCTIONS (FIVE PER FAULT CATEGORY). IT THEN BUFFERS IN 
C FROM FILE 'DBLFL' THE TWENTY-FIVE POSSIBLE DOUBLE-FAULT FUNCTIONS 
C (ALL COMBINATIONS OF PAIRED FAULT CATEGORIES). 

C THE SINGLE-FAULT FUNCTIONS ARE BUFFERED IN AND PLOTTED PER FAULT 
C CATEGORY. THE DOUBLE-FAULT FUNCTIONS ARE BUFFERED IN AND PLOTTED 
C PER FAULT CATEGORY PAIR. 

C NOTE - TAPE4=PLFILE IS USED BY THE DISSPLA ROUTINES. 

C 


COMMON/SNGFUNC/ PBNG ( 51 3 ) » PBGSTEP » NPBGSTP ( 20 ) » 

1 PNBNG(513) f PNBSTEP»NPNBSTP(20) r 

2 PFLD(513) »PFSTEP»NPFSTP(20) » 

3 PLAT(513) »PLTSTEP»NPLTSTP(20) » 

4 PDP(513) »PDPSTEP»NPDPSTP(20) f LNORLG 

C FOR OPTIMIZATION PURPOSES DURING THE BUFFER IN OF /SNGFUNC/r 
C DEFINE 'PLTSNG' TO ENCOMPASS THE ENTIRE COMMON AREA. 
DIMENSION PLTSNG(2671) 

EQUIVALENCE ( PBNG ( 1 ) » PLTSNG ( 1 ) ) 


EXTERNALS 

BUFBLK 

CPLOT 


TYPE ARGS 
6 
8 


DONEPL 0 
GENTMAR 4 
STPINDX 5 
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3.2.2 COVRGE and CVGPLT Files 



File Name 

Unit Name 


Description 

COVRGE Files; 




1. 

COVIN=^ , 

TAPE3 = 

COVIN, 

coverage definition input file 
generated by CAREIN program. 

2. 

INPUT , 

TAPES = 

INPUT 


3. 

OUTPUT , 

TAPE 6 = 

OUTPUT 


4 ! 

CVGMTS=J3 , 

TAPE? = 

CVGMTS, 

binary coverage moment functions 
passed to CARE3 program. 

5. 

DBUG, 

TAPES = 

DBUG 


e! 

SNGFL=jar, 

TAPES = 

SNGFL, 

plot file of single-fault prob- 
ability functions, passed to 
CVGPLT program. 

?! 

DBLFL=pr , 

TAPEIO = 

= DBLFL, 

plot file of double-fault prob- 
ability functions, passed to 
CVGPLT program. 

CVGPLT Files: 




1. 

COVIN=pf , 

TAPE 3 = 

COVIN , 

coverage input file required 
in CVGPLT also. 

* 

2. 

PLFILE , 

TAPE4 = 

PLFILE, 

DISSPLA plot file. 

3. 

INPUT, 

TAPES = 

INPUT 


4. 

OUTPUT , 

TAPE6 = 

OUTPUT 


5. 

SNGFL=j3f, 

TAPES = 

SNGFL, 

single-fault binary file from 
COVRGE program. 

6 . 

DBLFL=j? , 

TAPEIO = 

= DBLFL, 

double-fault binary file from 
COVRGE program. 


i* 

Must be saved and passed to other programs. 

*Must be saved and passed to DISSPLA plot program. 
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jzr 

Bannif/' 





C; f 

'loierc 

~D)OUi 

v/r 



CcmpLhk, 
Detnsity ^ 

fu/ichot) 's 
■&ro Foin4s 
Based ofl'W.u. 




SNGFLr 

fvL', -6 Ti/ichi^s 


pu4^a p(°r 
tsul-/- Ti/P£ 




■u/Llr Ti,Lt^ Tn/iic^ti/}s. 


(l£,yripu4t'C^p> £> 
^// Ty/’d 
f^l r 


3roF 


3.2.3 COVRGE Functional Flow Chart 


Jyrhoi Ccyn/r)CfJ 
/is/ fi^tofT)/ ^or 
C/9£o /YlcJu/e 


RAYTHEON COMPANY 
LtKINCTOM MASS 071?) 


"cmTg'L 


9956 














3.2.4 COVRGE Subroutines 


C 

C 

C 

C 


C 

C 


SUBROUTINE DBLFLT 


THIS SUBROUTINE COMPUTES THE DOUBLE-FAULT COVERAGE FUNCTION FOR 
COVERAGE FAULT TYPE PAIR ( ' I TYP ' r ' JTYP ' > * 


COMHON/WRKING/ B 1 AR ( 51 3 ) » BISTEP » NBl STP ( 20 ) » B2AR ( 51 3 ) » B2STEP » 
1 NB2STP(20) fDUM<513) rNDUM(20> »STPINIT(9) 

COMMON/FUNCCOM/ Cl AR ( 5 1 3 ) r CISTEP » NC 1 STP ( 20 ) , 

1 C2AR(513) f C2STEPrNC2STP(20) r 


2 

3 

4 

5 

6 

*********** STORAGE 


F1AR(513) tFISTEPjNFISTP 
F2AR(513) »F2STEP»NF2STP 
C3AR(513) »C3STEP»NC3STP 
C4AR(513) fC4STEP»NC4STP 
P3AR<513) fP3STEPrNP3STP 
AREA FOR FINAL DOUBLE- 


( 20 ) » 
( 20 ) » 
( 20 ) » 
( 20 ) » 
( 20 ) » 
FAULT 


FUNCTION ************ 


7 PDFAR(513) f PDFSTEPf NPDFSTP(20 ) »LNORLGj 
^ 1 / 

^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p 

8 XB1INT6(513)»XB1STEP»NXB1STP(20) » 

9 XB2INTG(513) fXB2STEP»NXB2STP(20) 


EXTERNALS TYPE ARCS 


BUFBLK 

6 

COMPFUN 

6 

FBDBL 

0 

FCDBL 

0 

FFDBL 

0 

GENMNTS 

4 

PREVNRC 

13 

SUMARS 

11 

TMAXDBL 

7 

VOLTERA 

10 


See the next p&9^ for the Double-fault Model Equations 
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DOUBLE-FAULT MODEL EQUATIONS 


aj (t) = ax e 
2 2 


-ait 

ax (t) = e ^ 


-ait 
2 


di (t) 
2 


rx(t) 

2 


= 1 


^1 (t) 
2 


bx (t) 
2 


-Bit 
Bx e ^ 


FUNCTION 

CNSRTDN 


■S 

-s 


6 1 (t) dT 
0 2 


p 1 (t) dx 

0 2 J 


FUNCTION 

RTDNINT 


= e 


-Bit 
2 


VARIABLE 


Cl (t) = Bi (t) da (t) ra (t) ag (t) + (1-P a)bi (t) 6a (t) ra (t) aa (t) 

2 2 111 ■/'I 


^1 2 


+ bi (t)da (t)Pa (t)aa (t) 
2 111 


CIAR(IT) 

C2AR(IT) 


f 2 (t) = ai (t)ba (t) di (t) ri (t) 

1 2 12 2 


C4(t) = 


^ ICi(t- 


T)Ba(T)bi(T) + Ca (t-T) Bi (T)ba (t) ]dx 


F2AR(IT) 

FIAR(IT) 


C4AR(IT) 


C 3 (t) 


■ L' 


x) Ba (T)bi (x) + f a (t-x) Bi (x)ba ( t) 1 dx 


C3AR(IT) 


s: 


P3 (t) = fi(t) +\ C3 (t-X)p3 (x)dx 


P3AR(IT) 


Df 


1 : 


(t) = Cl (t) +\ C4 (t-x)p 3 (x) dx 


PDFAR(IT) 
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SUBROUTINE COHPFUN < STEPS! » EXTFUN » ITH f T^M AX r NSTEPAR f FUNCAR ) 

C 

COMMON/CVGSTEP/ ITSTPS » MAXSTP » RELSTEP » TBASE » NTYPS > DBLDF » TRUHC » 
1 INMAXf NSTPMXf DIFCHN13 

COMhON/FLTYPCM/ I T YP » JTYP r CUPRNT 5 CVPLOT » I AXSTYP 
COMMON/CHICOM/ MDP » HLAT » HFLD f MBNG » MNBNG > MDF t CH I DBG 
DIMENSION FUNCAR(l) rHSTEPAR(l) 

LOGICAL THRASH fCVPRNTfCVPLOTfCHIDBG 

DATA DIVDIF/0* 1/f MAXDIV/4/»REALMIN/l *0E-293/ 

C 

C THIS SUBROUTINE COMPUTES A FUNCTION USING INCREASING STEP SIZESf 
C WHERE THE STEP SIZE IS DOUBLED EACH TIME IT INCREASES. THE 
C NUMBER OF STEPS COMPUTED FOR EACH STEP SIZE IS STORED IN NSTEPAR 
C ARRAY. 

C 


EXTERNALS 


TYPE ARCS 


EXTFUN REAL 2 F.P. 


INLINE FUNCTIONS TYPE 
ABS REAL 

AMAXl REAL 


ARGS 

1 INTRIN 
0 INTRIN 
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SUBROUTINE CULTAR ( SF AR » SFSTPST , NSFSTP ) 


C 

C COMPUTE UOLTERRA EQUATION USING SINGLE-FAULT ARRAY 'SFAR'» ITS 
C INITIAL STEP SIZE 'SFSTPST'f AND ITS 'NSFSTP' ARRAY DEFINING HOW 
C MANY OF EACH 'POWER OF TWO* TIMES 'SFSTPST' STEPS WERE GENERATED* 

C THESE ARRAYS WILL BE USED TO STORE THE FINAL RESULT OF THE CALL 
C TO SUBROUTINE 'VOLTERA' UPON RETURN. ' FAR ' » ' FSTPST ' AND 'NFSTPAR' 

C WILL BE USED AS TEMPORARY STORAGE FOR USE IN THE VOLTERA SUBROUTINE. 
C 

C 'FAR(TIME)' = 'SFAR(TIME)' + INTEGRATION FROM 0.0 TO TIME OF 
C 'FEEAR(TIME-TAU) ' * 'FAR(TAU)'. 

C 

COMMON/CVGSTEP/ I TSTPS f MAXSTP f RELSTEP r TBASE » NT YPS » DBLDF » TRUNC » 

1 INMAXtNSTPMXfDIFCHNG 

COMMON/WRKING/ FAR < 51 3 ) » FSTPST » NFSTPAR ( 20 ) , FEEAR ( 51 3 ) f FEESTEP f 
1 NFEESTP(20) »GAR(513) »NGSTPAR(20) fGSTPST(9) 

DIMENSION SFAR(513) »NSFSTP<20) 


EXTERNALS 

UOLTERA 
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o o n 


SUBROUTINE FILLDBL 


THIS SUBROUTINE FILLS THE DOUBLE-FAULT FUNCTIONS IN COHMON/CVRGCOH/ ♦ 


C 


COMHON/WRKING/ DFMOAR ( 513 ) » DFMSTEP » NDFMSTP ( 20 ) » DFHl AR ( 513 ) » 

1 DFM1DUM»NDFM1DH(20) »DFM2AR(513) »NDFM2DM(20) »DFM2DUM(9) 

COMMON/FLTYPCM/ I TYP » JT YP » CVPRNT » CVPLOT » I AXSTYP 
COMHON/CVGSTEP/ I TSTPS » MAXSTP » RELSTEP » TBASE f NT YPS » DBLDF f TRUNC > 


1 

1 

2 


INMAXrNSTPMXfDIFCHNG 
COMMON/CVRGCOM/ CMSTO ( 5 » 65 » 5 ) f CMST l ( 

CHDF0(5f5f65) fCMDFI ( 
TZER0ST(5f5) fT0DF(5f 
LOGICAL TRANSf CVPRNTfCVPLOT 
COMHON/TMARCOM/ DFTMAR ( 5 1 3 ) f TMAR2 ( 51 


5f65f5> fCMST2( 
5f5f65) fCHDF2( 
5) fTRANS(5) 

3) fTMAR3(513) f 


5f65f5) F 
5f5f65) F 


TMAR4<513) 


EXTERNALS TYPE 

FLINTP REAL 
GENTMAR 
STPINDX 


ARGS 

5 

A 

5 
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n o o 


SUBROUTINE FILLSNG ( MCHI ) 


THIS SUBROUTINE FILLS THE SINGLE-FAULT FUNCTIONS IN COMMON/CVRGCOM/ ♦ 


C 


1 


1 

1 


COMMON/WRKING/ SFMOAR ( 513 ) , SFMSTEP » NSFMSTP ( 20 ) » SFMl AR ( 513 ) » 

SFM1DUH»NSFH1DM(20) rSFH2AR(513) »NSFM2DM<20)fSFH2DUM(9) 
COMMON/FLTYPCM/ ITYP » JTYP » CVPRNT » CVPLOT i I AXSTYP 
COMMON/CVGSTEP/ ITSTPS » HAXSTP » RELSTEP » TBASE y NTYPS » DBLDF > TRUNC f 

INMAXfNSTPMXfDIFCHNG 

COMMON/CVRGCOM/ CMSTO ( 5 » 65 » 5 ) r CMST 1 ( 5 » 65 » 5 ) » CHST2 ( 5 f 65 r 5) » 

CHDF0(5»5r65) »CMDF1 (5»5»65) »CHDF2(5»5»65) » 
TZER0ST(5»5) »T0DF(5»5) »TRANS<5) 

LOGICAL TRANS, CMPRNTfCVPLOT 

COMMON/TMARCOM/ SFTMAR(513) , THAR2 ( 513 ) , TMAR3 ( 51 3) , TMAR4 ( 513 ) 


EXTERNALS TYPE ARGS 

FLINTP REAL 5 

GENTMAR 4 

STPINDX 5 
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n n n o 


SUBROUTINE GENMNTS ( FLT AR , FLTSTEP » NFLTSTP f HCHI ) 

C 

COMMON/CUGSTEP/ I TSTPS » MAXSTP . RELSTEP f TEASE » NTYPS » DBLDF » TRUNC r 
1 INriAX»NSTPMXfDIFCHNG 

COHHON/FLTYPCM/ I TYP » JTYP » CVPRNT » CVPLOT , I AXSTYP 
COMMON/WRKING/ FMOAR ( 513 ) » FHSTEP » NFHSTP ( 20 ) » FHl AR ( 5 1 3 ) » FM 1 BUM » 

1 NFM1DH<20) »FM2AR(513) »NFM2DM<20) »FM2DUH(9) 

LOGICAL CVPRNTrCVPLOT 
DIMENSION FLTAR(l) »NFLTSTP(1) 

DATA CNSONE/1.0/ 

IF THE FAULT FUNCTION IS CONSTANTLY ZERO» THERE IS NO NEED TO COM- 
PUTE THE MOMENT FUNCTIONS» SINCE THE MOMENT ARRAYS ARE INITIALIZED 
TO ZERO. 


EXTERNALS TYPE ARGS 

FILLDBL 0 

FILLSNG 1 

USTPINT 6 
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SUBROUTINE GENTMAR < STEPS! t NSTEPAR » THAR r TMX > 

C 

DIMENSION NSTEPAR(l) »TMAR(1) 

COMMON/CVGSTEP/ ITSTPS » MAXSTP » RELSTEP » TBASE » NTYPS » DBLDF r TRUNC » 
1 INMAX»NSTPMX»DIFCHNG 

COMHON/CHICOM/ MDP , HLAT » MFLD » MBNG » MNBND » MDF » CHIDBG 
LOGICAL CHIDB6 
C 

C THIS SUBROUTINE GENERATES AN ARRAY OF TIMES FOR WHICH A FUNCTION 
C WAS PREVIOUSLY COMPUTEIU 'STEPS!' IS THE INITIAL STEP SIZE USED 
C TO COMPUTE THE FUNCTION. 'NSTEPAR' IS THE ARRAY DEFINING HOW 
C MANY OF EACH ’POWER OF 2’ TIMES 'STEPS!' WERE USED. 

C I.E. NSTEPAR(NSTPIN) = THE NUMBER OF STEPS GENERATED HAVING 
C A STEP SIZE OF ’ 2** < NSTPIN-1 ) * STEPS!’. 

C 


59 



SUBROUTINE PREVNRC ( F » FSTEP f NFSTEP » FEE » FEESTEPi NFEEBTP , 

1 P»PSTEP»NPSTEP»CNSINTGf PAR»PARSTEP»NPARSTP) 

C 

COMHON/TMARCOM/ FTMAR(513) » FEETHAR ( 5 13 ) » PTMAR ( 51 3 ) » PARTHAR ( 513 > 

C 

C THIS SUBROUTINE'S MAIN FUNCTION IS AS A PREPROCESSOR FOR CALLING 
C SUBROUTINE 'VLTNREC'. THE TWO ARRAYS INVOLUED WITHIN THE CONVOLUTED 

C integration: 'FEE' and 'P' must be tested for their time maximums 

C (MAXIMUM X VALUE) AND PASSED TO SUBROUTINE 'VLTNREC' SO THAT THE 
C FUNCTION WITH THE LARGEST TIME MAXIMUM IS PASSED AS THE FUNCTION 
C THAT IS INTEGRATED FORWARD IN TIME* 

c example: 

C 'PAR'(T) = INTEGRATION FROM 0,0 TO 'T' OF 'FEE'(T-TAU) * 'P'(TAU), 

C IF 'P' FUNCTION HAS THE LARGEST TIME MAXIMUM (X VALUE) » THE FUNCTIONS 
C SHOULD BE PASSED IN THE SAME ORDER AS ABOVE. IF 'FEE' FUNCTION HAS 
C THE LARGEST TIME MAXIMUM^ THE EQUATION SHOULD BE AS FOLLOWS: 

C 'PAR'(T) = INTEGRATION FROM 0.0 TO 'T' OF 'P'(T-TAU) * 'FEE'(TAU). 


EXTERNALS 


GENTMAR 

VLTNREC 


TYPE ARGS 
4 

18 
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n n n 


SUBROUTINE PRNTCUG 
C 

ITSTPSf MAXSTP»RELSTEP»TBASE»NTYPS»nBLDF»TRUNC» 
INMAX fNSTPMX,DIFCHNG 

CMST0(5»65f 5) »CMST1 (5,65»5) » CMST2 ( 5 f 65 » 5 ) » 
CMDF0(5»5»65) fCHDFl (5»5»65> » CMDF2 ( 5 » 5 » 65 ) » 
TZER0ST(5»5) ,T0DF(5»5) »TRANS<5) 


PRINT CVRGCOM COMMON AREA. 


COMMON/CVGSTEP/ 

1 

C0MM0N/CVR6C0M/ 

1 

n 

LOGICAL TRANS 
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SUBROUTINE SN6FLT 


C 

C THIS SUBROUTINE COMPUTES THE SINGLE-FAULT COVERAGE FUNCTIONS FOR 
C COVERAGE FAULT TYPE 'ITYP'. 

C 

COMMON/URKING/ FAR ( 5 13 > » FSTPST » NFSTPAR ( 20 ) » FEEAR < 51 3 ) »FEESTEPr 
1 NFEESTP(20) >GAR(513 ) fNGSTPAR<20) fGSTPST(9) 

COMMON/FUNCCOM/ PA ( 513 ) » PASTEP » NPASTP ( 20 ) f 

1 PBl (513 ) fPB1STEPfNPB1STP<20) F 

2 PERR(513) fPERSTEPfNPERSTP(20) f 

C **)K**!|c***!|c FINAL AREA FOR STORAGE OF SINGLE-FAULT FUNCTIONS ********* 

3 PEAR(513) fPESTEPfNPESTP(20) f 

4 PNEAR(513) fPNESTEPfNPNESTP(20) F 

5 PFLD(513) fPFSTEPfNPFSTP(20) F 

6 PSIA(513) fPSASTEPfNPSASTP(20) F 

7 PDP(513) fPDPSTEPfNPDPSTP(20> fLNORLGf 

4^ 4^ 4^ 4^ 4^ 4^ 4^ 4^C 4^^ 4^ 4^ 4^ 4^ 4^ 4^ 4^ 4^C 4^ 4^ 4^ 4^ 4^ 4^ 

8 PSIB(513) fPSBSTEPfNPSBSTP(20) F 

9 PB2(513) fPB2STEPfNPB2STP(20) 




SINGLE-FAULT MODEL EQUATIONS 


4)(t) = ae 






XB<t) 


-6t 


Jo 


r ( t) d(T) dx 


P^(t) 

Ci 


Pj,(t) 


p^(t) 

e 


Pe(t) 


Pe(t) 


Pf(t) = 


-at 


= e r(t)d(t) + B I ij)(t-T)P (t)dT 


J 0 


r 

J 0 


<Ji(t) + 3 \ (j) (t-T)P, (T)dx 

0 


■s: 


-ax 


e p (x)d(x)e(t-x)dx + 3 \ 4)(t-x)P (x)dx 


r 

J 0 


-at 


= e p(t)d(t) + 3 \ (})(t-x)p (x)dx 


s; 


S' 


t 

e 6(t)r(t) + 3 \ (|) (t-x) p-(x) dx 

0 ® 


L-C, ( ‘ 

J 0 


(1-C) \ p (x) e (t-x)dx 
0 ® 


= C^ ?e (t) e (t-x)^-^-^ 


ae 


-(a+e) (t-x) 


a+ e 


) 


dx + P~(t) 




s 


'i'g(x)e"^^^”''^dx 




'^(x)dx + Pg 'l'g(x)dx 


VARIABLE 

FEEAR(IT) 

PA(IT) 

PBl (IT) 

PERR(IT) 

PEAR(IT) 

PNEAR(IT) 

PFLD(IT) 

PSIA(IT) 

PSIB(IT) 

PB2(IT) 

PDP(IT) 
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SINGLE FAULT MODEL EQUATIONS (CONT.) 


LET 


WITH 


VARIABLE 


Fx(t) = F^(t) + 




+ 


(1-Pg) BXg(t-x) ]Fx(t) dx 


F^(t) 


P^Ct) 



Pj,(t) 

Pp(t) 

F^(t) 


P^(t) 

a j 

Pj,(t) 

Xb 

p^(t) 

Pf(t) 


PgCt) = P^(t) + Pg(t) 

p„(t) + p-(t) PERMANENT FAULTS 

PL<t) = ® ® 

PgCt) TRANSIENT FAULTS 


PBNG(IT) 

PNBNG(IT) 

PLAT(IT) 
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SUBROUTINE STPINDX ( T ? STEPS! > NSTEPAR , STPSZE 7 INBX ) 

C 

COMHON/CUGSTEP/ ITSTPS > MAXSTP » RELSTEP » TBASE » NTYPS » DBLDF » TRUNC » 
1 INMAXf NSTPHXfDIFCHNG 

DIMENSION NSTEPAR(l) 

C 

C THIS SUBROUTINE COMPUTES THE STEP SIZE AND INDEX GIVEN A TIME 
C FOR A FUNCTION WITH ORIGINAL STEP SIZE OF 'STEPS!'. 'NSTEPAR' 

C DEFINES HOW MANY OF EACH 'POWER OF TWO' TIMES 'STEPS!' STEPS 
C WERE COMPUTED FOR THE FUNCTION. 

C 
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SUBROUTINE SUHARS ( ARl f STEPl ST » NSTPARl » CNl » AR2 f STEP2ST f NSTPAR2 » CN2 » 
1 ARSUM»STEPST»NSTEPAR> 

C 

C THIS SUBROUTINE SUMS TWO ARRAYS 'ARl' AND 'AR2' YIELDING 'ARSUM'. 

C 'STEPIST' AND 'STEP2ST' ARE ARRAYS' 'ARl' AND 'AR2' INITIAL STEP 
C SIZES RESPECTIVELY. 'NSTEFAR'» WHICH WILL CONTAIN THE NUMBER OF 
C STEPS OF EACH 'POWER OF TWO* TIMES 'STEPST' FOR 'ARSUM'» IS 
C CREATED IN THIS SUBROUTINE ALSO. 

C 'CNl' AND 'CN2' ARE CONSTANTS TO MULTIPLY EACH FUNCTION BY BEFORE 
C THE SUM. 

C 

DIMENSION ARl (1) » AR2 < 1 ) f ARSUM ( 1 ) » NSTPAR 1 ( 1 ) » NSTPAR2 ( 1) f NSTEPAR ( 1 ) 
COMMON/TMARCOM/ A 1 TMAR ( 51 3 ) » A2TMAR ( 5 1 3 ) . TMAR3 < 5 1 3 > , TMAR4 < 51 3 ) 
COMMON/CVGSTEP/ I TSTPS » MAXSTP r RELSTEP , TBASE » NT YPS , DBLDF » TRUNC f 
1 INMAX»NSTPMX»DIFCHNG 

COMMON/CHICOM/ MDP » MLAT » MFLD » MBNG » MNBNG » MDF » CHI DBG 
LOGICAL A1FNSH»A2FNSH» THRASH I CHI DBG 


EXTERNALS 

FLINTP 

GENTMAR 


TYPE ARCS 
REAL 5 

4 


STPINDX 


INLINE FUNCTIONS 
ABS 
AMAXl 
AMINl 


TYPE 

ARCS 

' 

REAX 

1 

TNTRIN 

REAL 

0 

INTRIN 

REAL 

0 

INTRIN 
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n n n n o 


SUBROUTINE TMAXDBL ( TCIMX » TC2MX ^ TFIMX > TF2MX » TBIMX » TB2MX » STPINIT ) 


COMPUTE TIME MAXIMUMS FOR THE SIX INITIAL DOUBLE-FAULT FUNCTIONSf 
TWO C FUNCTIONSr TWO F FUNCTIONS AND TWO B FUNCTIONSr BASED ON THE 
TIME MAXIMUMS FOR THE VARIOUS RATE FUNCTIONS, 


COMMON/FLTYPCM/ 

COMMON/RATECOM/ 

1 

2 

3 

LOGICAL CVPRNTr 
DIMENSION STPIN 


ITYP» JTYP»CVPRNT»CVPLOT» lAXSTYP 
ALPHA(5) ,BETA(5) »DELTA(5) »RH0<5) fEPSLN 
DELCDN(5) »RH0CDN<5) »EPSCDN(5) » 

PACNS<5) »PBCNS(5) »CCNS(5> » 

TALPMX(5) »TBETMX<5) f TDELMX < 5 ) » TRHOMX ( 5 
CVPLOT»DELCDN»RHOCDNfEPSCDN 
IT(9) 


(5) f 


) fTEPSMX 


5> 


INLINE FUNCTIONS 

TYPE 

ARGS 


AMAXl 

REAL 

0 

INTRIN 

AMINl 

REAL 

0 

INTRIN 
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noon 


SUBROUTINE THAXSNG ( TSN6MX » STPINIT ) 


COMPUTE TIME MAXIMUMS FOR THE NINE INITIAL SINGLE-FAULT G 
BASED ON THE TIME MAXIMUMS FOR THE VARIOUS RATE FUNCTIONS 


1 

o 

3 


COMMON/FLTYPCM/ ITYP » JTYP » CVPRNT » CVPLOT » lAXSTYP 
COMMON/RATECOM/ ALPHA<5) » BETA ( 5 ) » DELTA ( 5 ) »RH0(5> f EPSLN 

DELCDN(5) »RH0CDN(5) ,EPSCDN(5) , 

PACNS<5) ,PBCNS(5> »CCNS(5) » 

TALPMX(5) f TBETMX(5) ,TDELMX(5) >TRH0MX(5 
DIMENSION TSNGMX(9) ,STPINIT(9) 

LOGICAL CVPRNT » CVPLOT . DELCDN » RHOCDN » EPSCDN 


INLINE FUNCTIONS TYPE ARGS 

AMAXl REAL 0 INTRIN 

AMINl REAL 0 INTRIN 


FUNCTIONS 
5) T 

»TEPSMX(5) 
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SUBROUTINE VLTNREC ( F » FSTEP » NFSTEP » FTMAR » FEE » FEESTEP » NFEESTP f 

1 FEETMAR»P»PSTEpFNPSTEPfPTMAR»PTMX»CNSINTG»PARf PPRSTEP fNPARSTPf 

2 PARTHAR) 

C 

COMMON/CVGSTEP/ I TSTPS i HAXSTP r RELSTEP » TBASE » NTYPS » DBLDF » TRUNC f 
1 INMAXfNSTPMXfDIFCHNG 

COMHON/CHICOM/ MDP f MLAT f HFLUf MBNG f MNBNG f MDF f CHIDBG 
COMMON/TMPNTRS/ T f IFEEPNT f IPPNT f FEEFLG f PFLG 

DIMENSION F( 1 ) fFEE( 1 ) fP(1 ) fPAR( 1 ) fNFSTEP( 1 ) f NFEESTP < 1 ) f NPSTEP (1 ) f 
1 NPARSTP<1) fFTMAR(I) f FEETMAR ( 1 ) f PTHAR ( 1 ) f PARTMAR ( 1 ) 

LOGICAL THRASH fRECRSUFf FEEFLG F PFLG F CHIDBG 
DATA REALMIN/1 fOE-293/ 

C 

C THIS SUBROUTINE IS VERY SIMILAR TO THE UOLTERA SUBROUTINE* THE 
C ONLY DIFFERENCE IS THAT THIS SOLVES A NON-RECURSIVE CONVOLUTED 
C INTEGRAL WHERE ALL FUNCTIONS ARE KNOWN. 

C 'PAR'<T) = 'F'(T) + THE INTEGRATION FROM ZERO TO ' 1 ' OF 
C 'FEE'(T-TAU) * 'P'(TAU). 

C 


EXTERNALS 

TYPE 

ARGS 


CNVLINT 

REAL 

12 


FLINTP 

REAL 

5 


STPINDX 


5 


INLINE FUNCTIONS 

TYPE 

ARGS 


ABS 

REAL 

1 

INTRIN 

AMAXl 

REAL 

0 

INTRIN 

AMINl 

REAL 

0 

INTRIN 
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SUBROUTINE VOLTERA ( F » FSTEP » NFSTEP » FEE i FEESTEP f NFEESTP » CNS I NTG » 

1 P»PSTEPrNPSTEP> 

C 

COMMON/CVGSTEP/ I TSTPS » MAXSTP » RELSTEP » TBASE » NTYPS » DBLDF » TRUNC » 

1 INMAXfNSTPHX»DIFCHNG 

COMHON/CHICOM/ MDP » HLAT > MFLD» MBNG » MNBNG » MDF » CHIDBG 
COMMON/ThPNTRS/ T » I FEEPNT » I PPNT » FEEFLG » PFLG 
DIMENSION F(l) »FEE(1) »P(1) fNPSTEP(l) fNFEESTP(l) »NFSTEP(1) 
COMMON/TMARCOM/ FTMAR ( 5 1 3 ) r FEETMAR ( 51 3 ) » PTMAR ( 513 ) » TMAR4 ( 513 ) 
LOGICAL THRASH »STDYFLG»RECRSUF»FEEFLG,PFLG»CHIDBG 
DATA REALMIN/1 .OE-293/ 

C 

C UOLTERA SOLVES ANY LINEAR NON-SINGULAR VOLTERRA INTEGRAL EQUATION 
C OF THE SECOND KIND WITH CONVOLUTION KERNEL^ WHERE ' F ' f ' FEE '( KERNEL > 
C ARE GIVEN ARRAYS WITH INITIAL STEP SIZES 'FSTEP' AND 'FEESTEP', 

C 'F' IS THE COMPUTED ARRAY. THE LOWER LIMIT OF INTEGRATION IS 
C ASSUMED TO BE 0.0. 


EXTERNALS 

TYPE 

ARGS 


CNVLINT 

REAL 

12 


FLINTP 

REAL 

5 


GENTMAR 


4 


STPINDX 


5 


INLINE FUNCTIONS 

TYPE 

ARGS 


ABS 

REAL 

1 

INTRIN 

AMAXl 

REAL 

0 

INTRIN 

AMINl 

REAL 

0 

INTRIN 
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SUBROUTINE VSTPINT ( AR » STEPS! » NSTEPAR » CNSINTG » MOMNT , ARI NT ) 


C 

C THIS SUBROUTINE INTEGRATES A FUNCTION STORED IN 'AR' WHICH HAS A 
C GIVEN INITIAL STEP SIZE - 'STEPS!' AND A GIVEN SET OF STEPS WITH 
C INCREASING STEP SIZE. 'NSTEPAR' CONTAINS THE NUMBER OF STEPS GENER- 
C ATED FOR EACH 'POWER OF TWO* TIMES 'STEPS!'. 'ARINT' WILL HAVE THE 
C SAME 'STEPS!' AND 'NSTEPAR' AS 'AR'. 

C INTEGRATE THE ARRAY (Y VALUE) MULTIPLIED BY THE TIME (X VALUE) 


c 

RAISED TO 

THE ' 

MOMNT' POWER. 




c 

IF 

'MOMNT' 

= Of 

INTEGRATE 

THE 

ARRAY 

WITH NO MULTIPLIER. 

c 

IF 

'MOHNT' 

= 1» 

INTEGRATE 

THE 

ARRAY 

MULTIPLIED 

BY THE TIME. 

c 

IF 

'MOMNT' 

= 2» 

INTEGRATE 

THE 

ARRAY 

MULTIPLIED 

BY THE TIME SQUARED. 


C 

DIMENSION AR(1) rNSTEPARd ) f ARINT( 1) »REGINTG(20) 

COMMON/CVGSTEP/ ITSTPS » MAXSTP » RELSTEP » TBASE » NTYPS » DBLDF f TRUNC , 
1 INMAXfNSTPMXfDIFCHNG 

COMMON/CHICOM/ MDP » MLAT f MFLD » MBNG > MNBNG j MDF » CHI DBG 
LOGICAL CHIDBG 

COMMON/TMARCOM/ TMAR ( 513 ) r XFNATTM < 513 ) » TMAR3 ( 513) » TMAR4 ( 513 ) 

C 


EXTERNALS 

TYPE 

ARGS 

GENTMAR 


A 

SIMPINT 

REAL 

6 

STPINDX 


5 
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n n rj n o 


3.2.5 CQVRGE Functions 


FUNCTION ARMOMNK INDX»MOMNT»THARf AR) 

C 

niMENSION AR(1)»TMAR(1) 

THIS FUNCTION MULTIPLIES ' ( TMAR ( INDX ) **M0HNT ) *AR ( INDX > ' FOR INTEGRA- 
TION ROUTINE SIMPINT, 

THE TEST FOR A VALID MOMNT IS HANDLED IN SIMPINT. 
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FUNCTION CNSRTDN( RATE » TIME »DNFL6) 

LOGICAL DNFL6 
C 

C THIS FUNCTION RETURNS A CONSTANT RATE VALUE GIVE 'RATE' AND 'TIHE'» 
C IF DNFLG IS .FALSE, 


C 

C THIS FUNCTION RETURNS A CONSTANT DENSITY VALUE GIVEN 'RATE' AND 'TIME', 
C IF DNFLG IS .TRUE. 


EXTERNALS 

PREEXP 


TYPE ARCS 
REAL 1 
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FUNCTION CNVLINT(FVAL»FEE.FEETMAR»FEESTEPfNFEESTP, INDEXP»P»PTMAR» 
1 PSTEP»NPSTEP»RECRSVF»CNSINTG) 


C 

C THIS FUNCTION INTEGRATES FROM 0 TO 'T' (STORED IN COMMON/THPNTRS/ ) 
C CONVOLVING TWO FUNCTIONS STORED IN ARRAYS FEE AND P» 

C IF RECRSVFf RECURSIVE FLAG» IS .TRUEw THE P FUNCTION USED 
C IN THE INTEGRATION IS THE FUNCTION BEING COMPUTED. 

C 

COMMON/CVGSTEP/ I TSTPS i MAXSTP » RELSTEP , TBASE , NTYPS » DBLDF , TRUNC » 

1 INMAX»NSTPMX»DIFCHNG 

COMMON/TMPNTRS/ T f IFEEPNT » I PPNT » FEEFLG » PFLG 
LOGICAL PINTFLG»FEEINTF»FEEFLG»PFLG»RECRSVF 

DIMENSION FEEd ) »P(1) »FEETMAR(1) rPTMARd) »NFEESTP(1) »NPSTEP<1) 

C 


EXTERNALS 

TYPE 

ARGS 

FCNVLTM 

REAL 

2 

FLINTP 

REAL 

5 

FTCHSTP 

REAL 

3 

STPINDX 


5 
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FUNCTION FBDBL(ITHfT) 

C 

COMMON/RATECOM/ ALPHA ( 5 ) , BETA < 5 ) » DELTA ( 5 ) » RHO ( 5 ) , EFSLN ( 5 ) > 

1 DELCDN(5) »RH0CDN(5) »EPSCDN(5) » 

2 PACNS(5) »PBCNS(5> »CCNS(5) , 

3 TALPMX(5) ,TBETHX(5) »TDELMX(5) » TRHOMX ( 5 ) » TEPSMX ( 5 ) 
COMMON/FLTYPCM/ ITYP , JTYP » CUPRNT » CUPLOT » I AXSTYP 

LOGICAL DELCDN»RHOCDNfEPSCDN»DNSFLG»CVPRNT»CgpLOT 
DATA DNSFLG/. FALSE./ 

C 

C COMPUTE DOUBLE-FAULT B FUNCTION 1 OR 2. THE TWO FUNCTIONS^ IN- 
C VOLVING 'BETA'» USED TO DEFINE THIS FUNCTION ARE CONSTANT RATE 
C FUNCTIONS. 

C 


EXTERNALS TYPE ARGS 

CNSRTDN REAL 3 

RTDNINT REAL 3 
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c 


FUNCTION FCDBL(ITHfT) 


C 

C 

C 

C 

C 

C 


1 

2 

3 


COMMON/RATECOM/ ALPHA ( 5 ) , BETA( 5 ) » DELTA ( 5 ) » RHO ( 5 ) r EPSLN ( 5 ) » 

DELCDN(5) »RH0CDN(5> »EPSCDN(5) » 

PACNS(5) »PBCNS(5) f CCNS(5) » 

TALPMXO) fTBETHX(5) »TDELMX(5) »TRH0MX(5) »TEPSMX(5) 
COMMON/FLTYPCM/ I T YP » JT YP » CVPRNT f CVPLOT f I AXSTYP 
LOGICAL DELCDNfRHOCDNfEPSCDNfDNSFLGfCVPRNTfCVPLOT 


CONFUTE DOUBLE-FAULT C FUNCTION 1 OR 2, ALL FUNCTIONS USED TO COM- 
PUTE THIS C FUNCTION ARE CONSTANT RATE FUNCTIONS UNLESS FLAGS 
'RHOCDN' F 'DELCDN' » 'EPSCDN' SPECIFY CONSTANT DENSITY FOR FUNCTIONS 
INVOLVING 'RHO' F 'DELTA' OR 'EPSLN'* 


EXTERNALS TYPE 

CNSRTDN REAL 
RTDNINT REAL 


ARGS 

3 

3 
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FUNCTION FCNVLTM<FEETHAR»PTHAR) 

C 

DIMENSION FEETMAR(l) »PTMAR(1) 

COMMON/TMPNTRS/ T» IFEEPNT» IPPNTf FEEFLGf PFLG 
LOGICAL FEEFLGrPFLG 
C 

C FCNVLTM - FETCH CONVOLUTED TIME. 

C THIS FUNCTION RETURNS THE MAXIMUM TIME NOT YET FETCHED FROM EITHER 
C 'T-FEETMAR' OR PTMAR. IFEEPNT IS THE POINTER INTO FEETMARJ IPPNT IS 
C THE POINTER INTO PTMAR. FEEFLG AND PFLG ARE INTERPOLATION FLAGSf 
C IF BOTH ARE SET .TRUE.f THEY ARE USED AS FLAGS SPECIFYING THE 
C COMPLETION OF TIMES RETRIEVED FROM ONE OF THE ARRAYS. 

C 
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c 


FUNCTION FFDBK ITH»T) 


C 

C 

C 

C 

C 

C 


1 

2 

3 


CGHHON/RATECOH/ ALPHA ( 5 ) f BETA < 5 ) » DELT A ( 5 ) » RHO ( 5 ) » EPSLN ( 5 ) > 

DELCDN(5) »RH0CDN(5) »EPSCDN(5) » 

PACNS<5) f PBCNS(5) f CCNS(5) » 

TALPMX(5) »TBETHX(5) »TDELMX(5> fTRH0HX(5) f TEPSMX(5 
COHMON/FLTYPCH/ I TYP » JT YP » CVPRNT » CVPLOT , I AXSTYP 
LOGICAL DELCDNf RHOCDNjEPSCDN,DNSFLG»CVPRNTf CVPLOT 


) 


COMPUTE DOUBLE-FAULT F FUNCTION 1 OR 2. ALL FUNCTIONS USED TO 
COMPUTE THIS F FUNCTION ARE CONSTANT RATE FUNCTIONS UNLESS FLAGS 


'RHOCDN' » 'DELCDN' » 'EPSCDN' SPECIFY 
INVOLVING 'RHO' » 'DELTA' OR 'EPSLN', 


CONSTANT DENSITY FOR FUNCTIONS 


EXTERNALS TYPE 

CNSRTDN REAL 
RTDNINT REAL 


ARGS 

3 

3 
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on non 


FUNCTION FGSNGL( ITH»T) 


COMPUTE SINGLE-FAULT G FUNCTIONS 1 THROUGH 9. 

COMMON/FLTYPCM/ I TYP » JTYP » CVPRNT » CVPLOT » I AXSTYP 
COHHON/RATECOM/ ALPHA ( 5 ) » BETA < 5 > » DELTA ( 5 ) » RHO ( 5 ) » EPSLN ( 5 ) » 

1 DELCDN(5) ,RH0CDN(5) »EPSCDN(5) » 

2 PACNS<5) »PBCNS(5> »CCNS(5) » 

3 TALPMX(5) »TBETMX(5) »TDELMX(5) »TRH0MX(5) »TEPSMX(5) 
LOGICAL DELCDN»RHOCriN,EPSCDN»DNSFLG»CVPRNTfCUPLOT 
DATA DNSFLG/. FALSE./ 

ALL FUNCTIONS ARE CONSTANT RATE FUNCTIONS UNLESS FLAGS DELCDN f RHOCDN » 
C EPSCDN SPECIFY CONSTANT DENSITY FOR FUNCTIONS INVOLVING RHOiDELTA 
C OR EPSLN, 

C 


EXTERNALS TYPE 

CNSRTDN REAL 
PREEXP REAL 
RTDNINT REAL 


ARGS 

3 

1 

3 
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1 n n n 


FUNCTION FLINTP(T»STEP»ARTOINT, J,TIMEJ) 

C 

C0MM0N/CV6STEP/ ITSTPS » MAXSTP » RELSTEP » TBASE » NTYPS f DBLDF » TRUNC » 
1 INMAX>NSTPMX»DIFCHNG 

DIMENSION ARTOINT(l) 

THIS FUNCTION LINEARLY INTERPOLATES THE ARTOINT ARRAY WITH 
NON-UNIFORM STEP SIZES. 
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FUNCTION FTCHSTP(STEPSTf NSTEPARf INDX) 

C 

COMMON/CUGSTEP/ I TSTPS » MAXSTP » RELSTEP » TBASE , NTYPS » DBLDF > TRUNC » 
1 INMAXfNSTPHXf DIFCHNG 

DIMENSION NSTEPAR<1) 

C 

C GIVEN AN INDEX INTO AN ARRAY WITH NON-UNIFORM STEP SIZES 
C (NOTE: THE STEP SIZES ARE ’POWERS OF 2’ TIMES THE INITIAL 
C STEP SIZE - 'STEPST' ) fREFURN THE STEP SIZE ASSOCIATED WITH 
C THAT POINT. 

C NSTEF-AR CONTAINS THE NUMBER OF STEPS GENERATED FOR EACH 
C ’POWER OF 2’ TIMES 'STEPST'. 

C 
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n n 


FUNCTION RTDNINKRATEfTIMEf DNFL6) 

LOGICAL DNFLG 

THIS FUNCTION CONFUTES '1-THE INTEGRATION FROM 0 TO 'TIME' OF THE 
C CONSTANT RATE FUNCTION' EXPLICITLY! IF DNFLG IS .FALSE* 

C 

C THIS FUNCTION CONFUTES '1-THE INTEGRATION FROM 0 TO 'TIME' OF THE 
C CONSTANT DENSITY FUNCTION' EXPLICITLY! IF DNFLG IS .TRUE. 


EXTERNALS 

PREEXP 


TYPE ARCS 
REAL 1 
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FUNCTION SIMPINK ITFROMf ITTO»STPSZE»MOMNT»TMAR,AR> 


C 


C 

C 

C 

C 

C 


DIMENSION AR(1)»TMAR(1) 

COMHON/CUGSTEP/ I TSTPS » MAXSTP » RELSTEP , TBASE » NTYPS f DBLDF » TRUNC » 
1 INHAX»NSTPHX»DIFCHNG 


THIS FUNCTION 
(TIME**MOMNT) 
'STPSZE' OVER 


INTEGRATES THE FUNCTION 
FROM 'ITFROM' TO 'ITTO' 
THIS REGION. 


IN 'AR' MULTIPLIED BY 
WITH UNIFORM STEP SIZE 


EXTERNALS TYPE ARGS 

ARMOMNT REAL 4 

INLINE FUNCTIONS TYPE ARGS 

mod INTEGER 2 INTRIN 
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3.3 CARE 3 Module 


The CARE3 module consists of two main programs, CARE3 and 

RELPLT, and several subroutines and functions. Program CARE3 

calculates the system reliability by calculating the failed state 

* 

probabilities Q£(t) . (See the CARE III Final Report for a des- 
cription of the mathematical model used in the reliability model 
implementation.) Program RELPLT creates the plot file of the 
summary reliability functions. 


*J. J. Stiffler and L. A. Bryant, CARE III Final Report , 
Phase II , Preliminary, Raytheon Company, May 1981. 
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3.3.1 CARE 3 and RELPLT Main Programs 


PROGRAM CARE3 (RELIN=0»TAPE4=RELINf INPUT > OUTPUT »TAPE5=INPUT» 

1 TAPE6=0UTPUT»CV6MTS = 0»TAPE7=CU6MTB»DBUG»TAPE8=DBUG» rAPE9 = 0f 

2 TAPE 1 0 = 0 f TAPE 11 = 0» TAPE12 = 0 »GHFNCS»T APE 13 = GHFNCS»BXYFL = 0» 

3 TAPE14=BXYFLf FT 15Ff TAPE 15 = FT15Ff PLTFL = 0» TAPE 16=PLTFL» TAPE 17 = 0) 


C 

C 

C ************************* *************************************** »#!*:**** 

C* 

C* 

C* 

C* 

C* 

C* 

C* 

C* 

C* 

C* 

C* 

C* 

C* 


* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 

c 


title: COMPUTER-AIDED RELIABILITY ESTIMATION, VERSION THREE 

acronym: care hi 

MODEL design: DR. J.J.STIFFLER 

SOFTWARE DESIGN AND IMPLEMENTATION: LYNN A. BRYANT 
date: march 23,19S1 
DATE edited: 5-03-81 
EDITED by: LYNN A. BRYANT 


PURPOSE: 


CARE III IS A GENERAL-PURPOSE RELIABILITY 
ESTIMATION TOOL FOR FAULT-TOLERANT AVIONICS 
SYSTEMS. 


EXTERNALS 

TYPE 

ARGS 

BUFBLK 


6 

EOF 

REAL 

1 

FNCK 

REAL 

n 

RLSBRN 


0 
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CARES continued 


COMMON// G0RHSF(5f 20»<S5) , HDFPTS ( 1 00 » 1 00 > » I TBLKSZ 

note: this blank common area is also used to hold the LOADERf 

IN ORDER TO KEEP THE LOAD FIELD LENGTH REQUIREMENT THE 
SAME SIZE AS THE EXECUTION FIELD LENGTH* 

COMMON/STEPCOM/ I TSTPS » MAXSTP » RELSTEP » TBASE » NSTGRN , KWT t PSTRNC » 

1 PRCODE»RLPLOT»IAXSRL»SYSMNT>TBCF 

DIMENSION RECK4) »REC2<6> 

EQUIVALENCE ( I TSTPS » REC 1 < 1 ) ) r <NSTGRN» REC2 ( 1 ) ) 

C 

COMMON/CONFIG/ NCONVEC ( 20 ) » MSRVVEC < 20 ) » NFLTCAT ( 20 ) » MFLT YPE ( 5 » 20 ) » 

1 NSTGS» IFSTG»NTYPSrTRNSFC<5f 20) rMAXSTGrMAXCATrMAXI YFfMXSBRN 
DIMENSION REC3(163) 

EQUIVALENCE < NCONVEC ( 1 ) » REC3 ( 1 ) ) 

C 

COMMON/RATES/ 0MGA<5f 20) »RLAM(5f20) 

DIMENSION REC4(200) 

EQUIVALENCE ( OMGA ( 1 » 1 ) » REC4 ( 1 ) ) 

C 

COMMON/CNFGVAR/ LFLTVEC ( 20 ) » QLTSUM < 65 ) » PSTSUM ( 65 ) r QPSTSM ( 65 ) » 

1 QSUMWT(14) fQNOEFCTf PNOEFCT»PSTCOM»STOPARMf SUMK(3) fDBGFLG 

C 'PLTREL' WILL ENCOMPASS THE THREE ARRAYS PASSED TO PROGRAM RELPLT 
C FOR RELIABILITY PLOTS. 

DIMENSION PLTREL(195) 

EQUIVALENCE (QLTSUM(l) »PLTREL(1 ) ) 

C 

COMMON/NONLDEP/ AX I AR ( 5 . 20 » 65 ) » RXAR ( 20 » 65 ) 

DIMENSION SRPSTF<65»20) 

C AFTER EACH SUBRUN» 'SRPSTF' IS COMPUTED AND BUFFERED OUT TO TAPE17. 

C 'AXIAR' IS NO LONGER NEEDED AT THAT POINTf UNTIL THE NEXT SUBRUN - 
C WHERE IT IS CREATED USING THAT SUBRUN DATA. 

EQUIVALENCE ( AXIAR ( 1 , 1 i 1 ) f SRPSTF ( 1 » 1 ) ) 

C 

COMMON/CHICOM/ MDP » MLAT » MFLD t MBNG » MNBNG » MDF » CHIDBG 
COMMON/CVRGCOM/ CMSTO ( 5 i 65 » 5 ) » CMSTl ( 5 , 65 , 5 ) » CMST2 ( 5 j 65 r 5 ) y 

1 CMDF0(5r5y65) yCMDFl (5y5y65) y CMDF 2 ( 5 y 5 y 65 ) y 

2 TZER0ST(5y5) yT0DF(5y5) y TRANS (5) 

C FOR OPTIMIZATION PURPOSES DURING THE BUFFER IN OF /CVRGCOM/y 
C DEFINE AND EQUIVALENCE 'CVRGAR' TO ENCOMPASS THE ENTIRE COMMON AREA. 
DIMENSION CVRGAR(9805) 

EQUIVALENCE ( CMSTO ( 1 y 1 y 1 ) y CVRGAR( 1 ) ) 

C 

DIMENSION PSTFAR(65y70) yMNTRMV(70) yPS1FPD(65) y PSTSMW ( 1 4 ) y TI TLE ( 80 ) 
DIMENSION NSTGAR(35) 

C COVERAGE ARRAYS ARE NO LONGER NEEDED WHEN 'PSTFAR' IS BUFFERED IN 
C FROM TAPE17. 

EQUIVALENCE ( CMSTO ( 1 y 1 y 1 ) y PSTFAR ( 1 y 1 ) ) 

DIMENSION PSTBUF(4550) 

EQUIVALENCE (PSTFAR( 1 y 1 ) yPSTBUF( 1 ) ) 

C 

LOGICAL TRANS yTRNSFCyDBGFLGyCHIDBGyRLPLOT 
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PROGRAM RELPLK INPUT, OUTPUT »TAPE5= INPUT »TAPE6=0UTPUT» PLF I LE» 
1 TAPE4=PLFILE,RELIN=0,TAPE2=RELIN,PLTFL=0,TAPE16=PLTFL) 

C 


C* * 

C* title: COMPUTER-AIDED RELIABILITY ESTIMATION, VERSION THREE * 

c* acronym: care hi * 

C* MODEL design: DR* J.J*STIFFLER * 

C* SOFTWARE DESIGN AND IMPLEMENTATION: LYNN A. BRYANT * 

C* date: march 23,19S1 * 

C* DATE edited: 4-30-81 * 

C* EDITED by: LYNN A. BRYANT * 

C!|c * 

C* purpose: CARE III IS A GENERAL-PURPOSE RELIABILITY * 

C* ESTIMATION TOOL FOR FAULT-TOLERANT AVIONICS * 

C* SYSTEMS. * 

C* * 


^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p ^p 

C 

C THIS PROGRAM PLOTS THE TOTAL SYSTEM UNRELIABILITY - 'QLTSUM'j THE 
C TOTAL PERFECT COVERAGE UNRELIABILITY - 'PSTSUM'J AND THE SUM OF 
C THE TWO - 'OPSTSM', USING FILE PLTFL GENERATED IN PROGRAM CARE3. 

C NOTE - TAPE4=PLFILE IS USED BY THE DISSPLA ROUTINES. 

C 

COMMON/PLTCOM/ QLTSUM ( 65 ) , PSTSUM < 65 ) ,QPSTSM(65) 

DIMENSION PLTREL(195) 

EQUIVALENCE ( QLTSUM ( 1 ), PLTREL ( 1 > > 

C 

COMMON/STEPCOM/ ITSTPS , MAXSTP , RELSTEP , TBASE , NSTGRN , KWT , PSTRNC , 

1 PRCODE,RLPLOT,IAXSRL,SYSMNT,TBCF 

DIMENSION RECK4) ,REC2<6) 

EQUIVALENCE ( ITSTPS , RECl < 1 ) ) ,( NSTGRN , REC2 ( 1 ) ) 

C 


EXTERNALS 

BUFBLK 

CPLOT 

DONEPL 


TYPE ARCS 
6 
8 
0 
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3.3.2 CARE3 and RELPLT Files 



File Name 

Unit Name 


Description 

CARE3 Files: 




1. 

RELIN=j2(, 

TAPE 4 = 

RELIN, 

system definition input file 
generated by CAREIN program. 

2. 

INPUT, 

TAPES = 

INPUT 


3. 

OUTPUT, 

TAPE 6 = 

OUTPUT 


4. 

CVGMTS=j2f, 

TAPE? = 

CVGMTS, 

binary coverage moments file 
from COVRGE program. 

5. 

DBUG, 

TAPE 8 = 

DBUG 


6. 


TAPE 9 = 

0r 

HLAT internal binary file. 

7. 


TAPEljaf 

= 0, 

GFLD internal binary file. 

8. 


TAPE 11 

= 0. 

GNBNG internal binary file. 

9. 


TAPE 12 

= 0, 

HDF internal binary file. 

10. 

GHFNCS , 

TAPE 13 

= GHFNCS, 

debug file, under variable 
CHIDBG control, contains G 
and H functions computed 
using CVGMTS file. 

11. 

BXYFL=^, 

TAPE 14 

= BXYFL, 

binary critical pair results 
file created in program CAREIN. 

12. 

FT15F, 

TAPE 15 

= FT15F, 

system minterm file from 
CAREIN program. 

13l 

PLTFL=j2f, 

TAPE16 

= PLTFL, 

plot file of Q SUM, p* SUM, 
and Q+P* SUM passed to pro- 
gram RELPLT . 

14. 


TAPE 17 

= 0, 

function used to compute per- 


feet coverage function (com- 
puted per subrun) . 


^Must be saved and passed to program RELPLT. 


88 



Unit Name 


Description 


File Name 
RELPLT Files: 


1. 

RELIN=j2f, 

TAPE 2 = 

RELIN, 

system definition input file. 

2. 

PLFILE, 

TAPE 4 = 

PLFILE , 

DISSPLA plot file. 

3. 

INPUT, 

TAPES = 

INPUT 


4. 

OUTPUT, 

TAPE6 = 

OUTPUT 


5. 

PLTFL=jJ, 

TAPE16 = 

= PLTFL, 

plot file of Q SUM, P* SUM, 
and Q+P* SUM. 
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3.3.3 CARE3 Functional Flow Chart 
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3.3.4 CARE3 Subroutines 


SUBROUTINE ABCST ( MCHI r I T » ISTG , ICAT » JSTG » JCAT » N f FUNC ? 

1 ACOEF»BCOEF»CCOEF) 

C 

COHMON/CONFIG/ NCONUEC < 20 ) » MSRVUEC ( 20 ) » NFLTCAT ( 20 ) » MFLTYPE < 5 » 20 ) » 
1 NSTGS» IFSTG>NTYPS»TRNSFC<5»20) rMAXSTG jMAXCATfMAXTYPtMXSBRN 
COMMON/STEPCOM/ ITSTPS > HAXSTP r RELSTEP r TBASE f NSTGRN » KWT » PSTRNC f 
1 PRCODE»RLPLOT» lAXSRLf SYSMNTf TBCF 

COMMON/CHI COM/ MDP » ML AT » MFLD r MBNG » MNBNG r MDF » CH I DBG 
LOGICAL TRNSFCf CHIDBG»SYSMNTrRLPLOT 
C 

C THIS SUBROUTINE COMPUTES THE A»B AND C COEFFICIENTS THAT ARE 
C MULTIPLIED BY THE CMST FUNCTION ARRAYSf GOTTEN FROM THE COVERAGE 
C MODELf TO YIELD THE FHSFST AND FHDFST FUNCTIONS. 

C 


EXTERNALS 

FUNC 


TYPE ARGS 
REAL 6 F.P. 


INLINE FUNCTIONS TYPE 
mod INTEGER 


ARGS 

2 INTRIN 
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SUBROUTINE ARZERO ( AR . N » ALLO ) 


THIS SUBROUTINE TESTS 'AR' FOR A NON-ZERO VALUE. IF ALL VALUES IN 
THE ARRAY EQUAL ZERO» 'ALLO' IS SET .TRUE. 

LOGICAL ALLO 
DIMENSION AR(1) 


93 



SUBROUTINE BUFF I N ( MCH I , I T » I TB » I UN I T ) 

C 

C THIS SUBROUTINE FILLS EITHER GORHSF(G OR H SINGLE FAULT FUNCTION 
C FOR 'ITBLKSZ' NUMBER OF TIME STEP RECORDS) OR HDFPTS(H DOUBLE FAULT 
C FUNCTION PER TIME STEP) ARRAY DEPENDING UPON MCHI. lUNIT SPECIFIES 
C WHICH LOCAL FILE TO READ FROM. EACH BLOCK OF RECORDS ON 'TAPE 
C (lUNIT)' WAS CREATED FOR 'ITBLKSZ' NUMBER OF TIME STEPS 
C FOR ALL FAULT CATEGORIES (OVER ALL STAGES). THE DOUBLE FAULT FILE 
C CONTAINS ONLY ONE TIME STEP RECORD PER BLOCK. 

C 

COMMON/CHI COM/ MDP f MLAT » MFLD r MBNG » MNBNG > MDF r CHI DBG 
COMMON/CONFIG/ NCONVEC ( 20 ) » MSRVUEC ( 20 ) » NFLTCAT < 20 ) , MFLT YPE ( 5 » 20 ) y 
1 NSTGSy IFST6»NTYPSyTRNSFC(5»20) yMAXSTGfMAXCATyMAXTYPyMXSBRN 
LOGICAL TRNSFCyCHIDBGrSYSMNT 

COMMON// G0RHSF(5y20y65) y HDFPTS ( 1 00 y 1 00 ) y ITBLKSZ 
COMMON/STEPCOM/ ITSTPS y MAXSTP y RELSTEP y 1 BASE y NSTGRN y KWT y PSTRNC y 
1 PRCODEyRLPLOTy lAXSRLySYSMNTyTBCF 

C 


EXTERNALS TYPE ARCS 

UNIT REAL 1 
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SUBROUTINE BUFFOUT ( MCHI » ITB f lUN 1 T ) 


C 

C THIS SUBROUTINE CREATES A BLOCK OF RECORDS ON LOCAL "^ILE 
C 'TAPE(IUNIT) ' USING BUFFER OUT» FROM EITHER GORHSF OR HDFPTS ARRAY 
C DEPENDING UPON MCHI* EACH BLOCK OF RECORDS CONTAINS 'ITBLKSE' 

C NUMBER Or TIME STEPS FOR ALL FAULT CATEGORIES <OVER ALL STAGES) 

C FROM GORHSF ARRAY. WHEN CREATING THE DOUBLE FAULT FILE USING 
C HDFPTS ARRAY» EACH BLOCK CONTAINS JUST ONE TIME STEP RECORD. 

C 

COMMON/CHICOM/ MDP > MLAT » MFLD » MBNG » MNBNG » MDF t CHI DBG 
COMMON/STEPCOM/ I TSTPS » MAXSTP f RELSTEP » TBASE » NSTGRN » KWT f PSTRNC » 

1 PRCODEfRLPLOTfIAXSRLjSYSMNTfTBCF 

COMMON/CONFIG/ NCONVEC ( 20 ) f MSRUVEC ( 20 ) f NFLTCAT ( 20 > f MFLT YPE ( 5 f 20 ) f 
1 NSTGSf IFSTGfNTYPSfTRNSFC(5f20) fMAXSTGfMAXCATfMAXTYPfMXSBRN 
LOGICAL TRNSFCfCHIDBGfSYSMNTfRLPLOT 
COMMON// G0RHSF(5f20f65) fHDFPTS<100f 100) fITBLKSZ 
C 


EXTERNALS 

UNIT 


TYPE ARGS 
REAL 1 
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SUBROUTINE BXYC( ISTGi JSTG» INBXY»LXC»LYC) 


C 

C FOR A 20 STAGE SUBRUN» 210 UNIQUE BXY FUNCTION DEFINITIONS ARE 
C POSSIBLEf BUT ONLY 55 DEFINITIONS CAN BE IN MEMORY AT ONE TIME, 

C 'IBREC' SPECIFIES WHICH RECORD IS CURRENTLY IN MEMORY. (A MAXIMUM 
C OF FOUR RECORDS CAN EXIST PER 20 STAGE BLOCK OF DEFINITIONS.) 

C 

C INDEX INTO 'BXYAR'5 BXYAR ( LX-MUX+1 » LY-MUY+1 » INBXY ) » WHERE LX AND lY 
C ARE THE CURRENT NUMBER OF FAULTS IN 'IST6' AND 'JST6' RESPECTIVELY 
C AND MUX=0.,LX AND MUY=O..LY. 

C 

C THIS SUBROUTINE COMPUTES 'INBXY' GIVEN THE 'IST6»JST6' PAIR. 

C 'IJSTGIN' CONTAINS THE INDEX 'INBXY' BASED ON THE CONVERSION OF 
C 'ISTGfJSTB' INTO 'IJSTG'. NOTE: THE BXY FUNCTION IS THE SAME FOR 
C THE PAIRS 'ISTGfJSTG' AND 'JSTGflSTG'r I.E. STAGE PAIR ORDER IS NOT 
C IMPORTANT. 

C IF 'INBXY' EQUALS ZERO» NO BXY FUNCTION EXISTS FOR THAT STAGE PA3K. 

C IF 'INBXY' IS GREATER THAN ' I BREC ' * ' MXBDEF ' » BUFFER IN THE NEXT 
C RECORD^ WHICH CONTAINS THE NEXT 55 DEFINITIONS OF THE BXY FUNCTIONS. 

C IF 'INBXY' IS .LE. TO ' IBREC ' * ' MXBDEF ' f BACKSPACE TWO RECORDS AND 
C BUFFER IN THE DATA. THIS IS POSSIBLE BECAUSE BXY FUNCTIONS ARE 
C USED IN BOTH ROUTINES 'FBCRTL' AND 'FDSCRTL' PER FAULT VECTOR. 

C note: in most CASES» only one record is required to DEFINE THE NON- 
C ZERO BXY FUNCTIONS FOR 'NSTGS' COUPLED STAGES SUBRUN. I.E. THERE 
C WILL ONLY EXIST 'MXBDEF' OR LESS NON-ZERO BXY FUNCTIONS. THEREFORE 
C ONLY THE INITIAL BUFFER IN WILL BE REQUIRED PER SUBRUN» WHICH IS DONE 
C IN SUBROUTINE NFLTVDP, 

C 

COMMON/CONFIG/ NCONVEC ( 20 ) » MSRVVEC ( 20 ) f NFLTCAT ( 20 > » MFLTYPE ( 5 » 20 ) » 

1 NSTGS»IFSTG»NTYPSf TRNSFC(5»20) »MAXSTG»MAXCAT,MAXTYP,MXSBRN 
LOGICAL TRNSFC»FATALfEOFFLG 

COMMON/BXYCOM/ BXYAR ( 1 0 » 10 » 55 ) » IJSTG IN ( 21 0 ) »LCNDVEC(20) f IBRECr 
1 KFSTG»KPLSTG(70) 

C FOR OPTIMIZATION PURPOSES DURING THE BUFFER IN OF /BXYCOM/» DEFINE 
C AND EQUIVALENCE 'BCOMAR' TO ENCOMPASS THE ENTIRE DATA TO BE BUFFERED. 
DIMENSION BC0MAR(5802> 

EQUIVALENCE ( BXYAR < 1 f 1 » 1 ) f BCOMAR ( 1 ) ) 


EXTERNALS 

TYPE 

ARGS 

BUFBLK 


6 

INLINE FUNCTIONS 

TYPE- 

ARCS — ^ 

MOD 

INTEGER 

2 INTRIN 
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SUBROUTINE CAXLAT 
C 

COMMON/CONFIG/ NCONVEC ( 20 ) » MBRVVEC ( 20 ) , NFLTCAT ( 20 ) , MFLTYPE ( 5 » 20 ) t 
1 NSTGS» IFSTG»NTYPS»TRNSFC(5»20) »MAXSTG»MAXCAT,MAXTYP»MXSBRN 
COHMON/CHICOH/ MDP » MLAT , MFLD f MBNG » MNBNG » MDF » CHI DBG 
COMMON/NONLDEP/ AXI AR ( 5 » 20 r 65 ) » RXAR ( 20 » 65 ) 

COMMON/RATES/ 0MGA(5»20) »RLAM(5»20) 

COMMON/STEPCOH/ I TSTPS » MAXSTP » RELSTEP > TBASE r NSTGRN r NUT f PSTRNC » 

1 PRCODEjRLPLOT» lAXSRLfSYSMNT f TBCF 

LOGICAL TRNSFC»CHIDBG»SYSHNT»RLPLOT 
COMMON// G0RHSF<5»20r65) »HDFPTS<100f 100) r ITBLKSZ 
DIMENSION HLAT<5f 20»65) 

EQUIVALENCE ( HLAT ( 1 » 1 » 1 ) » GORHSF ( 1 » 1 » 1 ) ) 

C 

C THIS SUBROUTINE COMPUTES ARRAY AXI AR ( I CAT » I STG » I T ) WHICH IS THE 
C PROBABILITY THAT A GIVEN STAGE X MODULE <ISTG) HAS A CATEGORY 
C XI (ICAT) LATENT FAULT AT TIME ’T'f REPRESENTED BY 'IT*, 

C GIVEN THAT IT HAS EXPERIENCED SOME FAULT BY TIME 'T'. 

C 


PERMANENT 

TRANSIENT 


f! xsra^ tt) 


EXTERNALS 

BUFFOUT 

FHSFST 

PRNTGH 


TYPE ARGS 
3 

REAL 4 

3 


c/£/n^ 

and 


FlJsFsr(?f)L/)T, 
Tsr^; .xcflr) 


INLINE FUNCTIONS TYPE 
MOD INTEGER 


ARGS 

2 INTRIN 
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SUBROUTINE CRXFF 
C 

COMMON/CONFIG/ NCONUEC ( 20 ) » MSRUUEC ( 20 ) »NFLTCAT(20) » MFLT YPE ( 5 » 20 ) » 

1 NSTGS»IFSTG»NTYPS»TRNSrC<5f 20) »MAXSTGf MAXCAT»MAXTYP»«XSBRN 
COMMON/NONLDEP/ AXI AR ( 5 » 20 » 65 ) » RXAR ( 20 » 65 ) 

COMMON/STEPCOM/ I TSTPS » MAXSTP > RELSTEP > TB ASE » NST6RN » KUT » PSTRNC > 

1 PRCODEjRLPLOT»IAXSRL>SYSHNTjTBCF 

COMMON/CHICOM/ MDP » HLAT » MFLD » MBNG » MNBNG » MDF » CHI DBG 
LOGICAL TRNSFC>CHIDBG»TRNSFLTfSYSMNT»RLPLOT 

THIS SUBROUTINE COMPUTES THE PROBABILITY THAT A SINGLE STAGE X MODULE 
IS FAULT-FREE AT TIME T» REPRESENTED BY ITf FOR ALL STAGES, 


RELIABILITY OF A STAGE x MODULE 




n r (t) 

i 


R.%f)R(rr) 


IT 

Xcftr 


TRXS.fLT') 


EXTERNALS 

FRXIFF 

PRNT6H 


TYPE ARCS 
REAL 4 

3 


98 



o n 


SUBROUTINE GNFLTVC ( I ST I NI T » ISTSTOP » MAXLAST ) 


C 

C 

C 

C 


THIS SUBROUTINE GENERATES FAULT VECTORS IN SETS STARTING AT 
SET ISTINIT AND STOPPING AT SET ISTSTOP. IF I ST I N I T . EQ . 1 r 
THE UNRELIABILITY OF THE SYSTEH IS BEING COMPUTED. IF ISTINIT 
.GT.lf PERFECT COVERAGE PROBABILITIES WILL BE COMPUTED. 


COMMON/CONFIG/ NCONVEC ( 20 ) » MSRVVEC ( 20 ) »NFLTCAT(20) , MFLT YPE ( 5 » 20 ) » 
1 NSTGS» IFST6»NTYPS»TRNSFC(5»20> »MAXST6»MAXCAT»MAXTYPf MXSBRN 


1 

1 

1 


COMMON/CNFGVAR/ LFLTVEC ( 20 ) » QLTSUM ( 65 ) » PSTSUM ( 65 ) » QPSTSM ( 65 ) 7 

QSUMWT( 14) 7QN0EFCT»PN0EFCT»PSTC0M»ST0PARM7SUMK(3) 7 DBGFLG 
COHMON/STEPCOM/ ITSTPS » MAXSTP 7 RELSTEP » TBASE 7 NSTGRN 7 KWT > PSTRNC r 

PRCODE 7 RLPLOT 7 IAXSRL 7 SYSMNT 7 TBCF 
COMMON/GNVCCOM/ QLT(65) 7 PSTLAR ( 65 ) 7 PLT ( 65 ) 7QINTGRL(3 ) 7 

ISETVC(20) 7lSETVCS(20) 7LVLAST(20) 7 QSUMSF 7 PSUMSF 
LOGICAL TRNSFC 7 DBGFLG 7 COMPQ 7 COMPPST 7 RLPLOT 
LOGICAL ONOEFCT7pNOEFCT7pSTCOMf ALLO 7 SYSMNT 


EXTERNALS 

TYPE 

ARGS 


ARZERO 


3 


FPSTAR 

REAL 

2 


UNRELQ 


2 


INLINE FUNCTIONS 

TYPE 

ARGS 


AMAXl 

REAL 

0 

INTRIN 

MINO 

INTEGER 

0 

INTRIN 
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SUBROUTINE NFLTVDP 


C 

C THIS SUBROUTINE AND SUBROUTINES CALLED WITHIN COMPUTE ALL FUNCTIONS 
C THAT ARE NOT DEPENDENT UPON THE FAULT VECTORS. 

C FUNCTIONS THAT ARE USED CONTINUALLY THROUGHOUT THE ROUTINES ARE STORED 
C IN COMHON/NONLDEP/ AXI AR ( 5 r 20 » ) » RXAR ( 20 » 65 ) . 

C IN ORDER TO SAVE MEMORYf CERTAIN G AND H FUNCTIONS^ NAMELY HLAKLATENT 
C FAULT) » GFLD(FAILED) » GNBNG(NOT BENIGN) AND HDF<DOUBLE FAULT) j ARE 
C COMPUTED IN THIS SUBROUTINE AND BUFFERED OUT TO A LOCAL FILE WHICH 
C WILL BE BUFFERED IN WHEN REQUIRED. LOCAL FILES TAPE9r TAPEIO » TAPEl 1 > 

C TAPE12 CONTAIN THESE G AND H FUNCTIONS. 

C TAPE13 NAMED GHFNCB IS A DEBUG PRINTOUT OF ALL G AND H FUNCTIONS 
C COMPUTED - CREATED IN SUBROUTINE PRNTGH IF CHIDBG IS .TRUE. 

C 

COMMON/CHICOM/ MDP » MLAT > MFLD » MBNG . MNBNG , MDF , CHI DBG 
COMMON/STEPCOM/ I TSTPS » MAXSTP » RELSTEP , TBASE » NSTGRN » KWT » PSTRNC s- 
1 PRCODE»RLPLOT»IAXSRL,SYSMNT»TBCF 

COMMON/NONLDEP/ AX I AR ( 5 » 20 f 65 ) » RXAR ( 20 , 65 ) 

COMMON/CONFIG/ NC0NVEC(20) » MSRVVEC ( 20 ) f NFLTCAT ( 20 ) » MFLTYPE < 5 » 20 ) , 

1 NSTGSf IFSTG»NTYPS»TRNSFC(5 t 20) » MAXSTG , MAXC AT r MAXT YP » MXSBRN 
LOGICAL TRNSFCt CHIDBG >SYSMNT,NONFTL»EOFFLG,RLPLOT 
COMMON// G0RHSF(5»20»65) »HDFPTS< 100» 100) » ITBLKSZ 
DIMENSION HBNG(5f 20»65) 

EQUIVALENCE (HBNG( 1 ,1,1 ) ,GORHSF( 1 ,1,1)) 

C FOR OPTIMIZATION PURPOSES DURING THE BUFFER IN OF /BXYCOM/, DEFINE 
C AND EQUIVALENCE 'BCOMAR' TO ENCOMPASS THE ENTIRE DATA TO BE BUFFERED. 
COMMON/BXYCOM/ BXYAR ( 10 , 10 , 55 ) , IJSTGIN ( 210 ) , LCNDVEC ( 20 ) ,IBREC, 

1 KFSTG,KPLSTG(70) 

DIMENSION BC0MAR(5802) 

EQUIVALENCE ( BXYAR ( 1 ,1,1) , BCOMAR < 1 ) ) 


EXTERNALS 

TYPE 

ARGS 

BUFBLK 


6 

BUFFOUT 


3 

CAXLAT 


0 

CRXFF 


0 

FGST 

REAL 

A 

FHDFST 

REAL 

6 

FHSFST 

REAL 

4 

FNCK 

REAL 

2 

PRNTGH 


3 

INLINE FUNCTIONS 

TYPE 

ARGS 

MOD 

INTEGER 

2 INTRIN 
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SUBROUTINE PRNTGH < HCHI » I UNI T t GORH ) 


C 

C THIS SUBROUTINE IS USED TO PRINT G AND H FUNCTIONS CURRENTLY STORED 
C IN COMMON//. MCHI DETERMINES WHETHER THE ARRAY TO BE PRINTED 
C IS A SINGLE OR DOUBLE FAULT ARRAY. WHEN MCHI.EQ.MDFt IUNIT SPECIFIES 
C WHICH FILE CONTAINS THE DOUBLE FAULT FUNCTION. SINCE THE DOUBLE 
C FAULT FUNCTION IS TOO LARGE TO STORE IN CENTRAL MEMORY» IT MUST BE 
C BUFFERED IN FROM DISK IN ORDER TO BE PRINTED. 

C 

COMMON// G0RHSF(5»20»65) » HDFPTS < 1 00 » 1 00 ) , I TBLKSZ 
COMMON/STEPCOM/ ITSTPS » MAXSTP » RELSTEP » TBASE » NSTGRN» KWT » PSTRNC r 
1 PRCODE»RLPLOTf IAXSRL>SYSMNT»TBCr 

COMMON/CHICOM/ MDP » MLAT » MFLD > MBNG » MNBNG » MDF » CHI DBG 
COMMON/CONFIG/ NCONVEC ( 20 ) » MSRVVEC ( 20 ) » NFLTCAT ( 20 ) » MFLT YPE ( 5 » 20 ) » 

1 NSTGS»IFSTG»NTYPS»TRNSFC(5»20) »MAXSTG»MAXCAT»MAXTYP»MXSBRN 
LOGICAL TRNSFC»CHIDBG»SYSMNT»RLPLOT 
DIMENSION GHSFPTS(5»20) 

EQUIVALENCE ( GHSFPTS < 1 ^ 1 ) f GORHSF ( 1 » 1 » 1 ) ) 

DIMENSION MCHICD<6> 

DATA MCHICD/4HDP »4HLAT MHFLD »4HBNG »4HNBNG»4HDF / 


EXTERNALS 

BUFFIN 


TYPE ARCS 
4 
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SUBROUTINE RLSBRN 


THIS SUBROUTINE COMPUTES THE SYSTEM UNRELIABILITY PER FAULT VECTOR 
FOR A SUBSET OF THE TOTAL NUMBER OF STAGES. A MAXIMUM OF 20 STAGES 
PER SUBRUN ARE USED. 


COMMON/STEPCOM/ I TSTPS » MAXSTP » RELSTEP t TBASE » NSTGRN » KWT » PSTRNC » 

1 PRCODEfRLPLOT» lAXSRLf SYSMNTf TBCF 

COMMON/CONFIG/ NCONVEC ( 20 ) » MSRVVEC ( 20 ) t NFLTCAT ( 20 ) » MFLTYPE < 5 » 20 ) » 
1 NSTGSf IFSTG»NTYPS.TRNSFC(5»20 ) jMAXSTGtMAXCATtMAXTYP^MXSBRN 
COMMON/RATES/ OMGA < 5 f 20 ) r RLAM ( 5 » 20 > 

COMMON/CNFGVAR/ LFLTVEC ( 20 ) r QLTSUM ( 65 ) f PSTSUM ( 65 ) » QPSTSM ( 65 ) » 

1 QSUMWT(14) TQNOEFCTrPNOEFCTjPSTCOMf ST0PARM>SUMK(3) fDBGFLG 

COMMON/NONLDEP/ AXI AR < 5 » 20 » 65 ) » RXAR ( 20 » 65 ) 

COMMON/CVRGCOM/ CMSTO ( 5 » 65 » 5 ) » CMST 1 ( 5 » 65 » 5 ) > CMST2 ( 5 » 65 » 5 ) , 

1 CMDF0<5»5»65) »CMDF1 (5»5»65) f CMDF2(5»5»65) » 

2 TZER0ST(5f5) »T0DF(5f 5) »TRANS(5) 

LOGICAL TRANS »TRNSFCjDBGFLG»RLPLOT 

LOGICAL QNOEFCT»PNOEFCT»PSTCOM»SYSMNT 


EXTERNALS 

TYPE 

ARCS 

GNFLTVC 


3 

NFLTVDP 


0 

INLINE FUNCTIONS 

TYPE 

ARGS 

MINO 

INTEGER 

0 INTRIN 
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SUBROUTINE SUMHAT < I T » NFLTS > 


C 

C 

C 

C 

C 

C 

C 

C 


COHMON/CONFIG/ NCONVEC ( 20 ) » HSRUVEC ( 20 ) » NFLTCAT ( 20 ) , HFLT YPE ( 5 » 20 ) , 
1 NSTGS»IFSTG»NTYPS»TRNSFC(5»20) »MAXSTGrMAXCATfMAXTYP»MXSBRN 
COMMQN/CNFGVAR/ LFLTVEC ( 20 ) t QLTSUM ( 65 ) » PSTSUM ( 65 ) » QPSTSH ( 65 ) » 

1 QSUMUT(14) fQN0EFCT»PN0EFCT»PSTC0M»ST0PARMfSUMK(3> »DBGFLG 

CDMMON/CHICOM/ MDP r MLAT f MFLD » MBNG r MNBNG f MDF r CHIDBG 
COMMON/GNOCCOH/ QLT ( 65 ) » PSTLAR ( 65 ) »PLT(65) ,0INTGRL(3) f 
1 ISETVC<20 ) fISETVCS(20) fLVLAST(20) t QSUMSF » PSUMSF 

COMMON/STEPCOH/ I TSTPS » MAXSTP f RELSTEP f TBASE f NSTGRN » KWT f PSTRNC » 


1 PRCODEfRLPLOTf IAXSRLfSYSMNTfTBCF 

COMMON// G0RHSF(5f20f65) fHDFPTS(IOOfIOO) F ITBLKSZ 
COMMON/BXYCOM/ BXYAR ( 1 0 f 1 0 f 55 ) f IJSTGI N ( 2 1 0 ) f LCNDVEC ( 20 ) f I BREC f 


1 KFSTGfKPLSTG(70) 

LOGICAL TRNSFCfDBGFLGfCHIDBGfQNOEFCTfPNOEFCTfPSTCOMfSYSMNTfRLPLOT 
DIMENSION LMUNTVC(20) 


THIS SUBROUTINE COMPUTES THE SUMMATIONf OVER ALL FAULT CATEGORIEHf 
OF ALL THE FUNCTIONS INVOLVED IN COMPUTING THE SUBSYSTEM RELIABILITY 
MODEL* THIS SUMMATION IS COMPUTED OVER ALL TIME STEPS FOR THE 
CURRENT LFLTVEC AND STORED IN SUNK ARRAY. SUMK IS INTEGRATED USING 
THE SIMPSON'S 1/3 RULE IN ORDER TO GENERATE THE OLT ARRAY FOR THE 
CURRENT LFLTVEC PER TIME STEP. 


K,(t) 


E 




+ A (t 1 

A)pJ(t) + 

a (t| 

i)P*(T) 

EXTERNALS 

TYPE 

ARCS 


FAPC 

REAL 

*m 


FCYJ 

REAL 

4 


FLAM 

REAL 

3 


FPSTAR 

REAL 

n 

*m 


FPSTREC 

REAL 

A 


INLINE FUNCTIONS 

TYPE 

ARGS 


MOD 

INTEGER 

n 

Am 

INTRIN 
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SUBROUTINE UNRELQ ( NFLTS » NFLSTG ) 

COHHON/GNVCCOM/ QLT ( 65 ) » PSTLAR ( 65 ) r PUT (65) »QINTGRL<3) f 
1 ISETVC(20) »ISE1UCS(20) >LVLAST(20) » QSUMSF » PSUHSF 

COMMDN/STEPCOM/ I TSTPS » MAXSTP » RELSTEP » TB ASE » NSTGRN » KWT , PSTRNC » 

1 PRCODE»RLPLOT»IAXSRL»SYSMNT»TBCF 

- COMHON/CNFGVAR/ LFLTVEC ( 20 > » QLTSUM ( 65 ) »PSTSUM(65) »QPSTSM(65) t 
1 QSUMWK 14) »QN0EFCT»PN0EFCT»PSTC0M»ST0PARM»SUMK(3) tDBGFLG 

LOGICAL QNOEFCT»PNOEFCT»PSTCOM»DBGFLG»SYSHNTf RLPLOT 

THIS SUBROUTINE COMPUTES THE SYSTEM UNRELIABILITY OVER ALL TIME 
STEPS PER FAULT VECTOR AS THEY ARE GENERATED IN SUBROUTINE 
GNFLTVC. 


MODIFIED RELIABILITY EQUATION 


Qj(t) = t Kjj(-T)dT 


Uaria Llf S 

QLT{C,S) 

^umkCa) 


EXTERNALS TYPE ARGS 

FINTGRT REAL 2 

SUMMAT o 
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3.3.5 CARE3 Functions 


FUNCTION FAC(IT»LVECTOR) 

C 

COMMON/CONFIG/ NC0NVEC<20) » «SR VVEC ( 20 ) r NFLTCAT ( 20 ) » MFLTYPE < 5 » 20 ) r 
1 NSTGSf IFST6fNTYPS»TRNSFC<5f 20) »MAXSTG»MAXCAT»MAXTYP»MXSBRN 
COMMON/CHICOM/ MDP » MLAT »MFLD »MBNG » MNBNG > MDF r CHIDBG 
COMMON/NONLDEP/ AXI AR ( 5 » 20 » 65 ) » RXAR ( 20 r 65) 

LOGICAL TRNSFCrCHIDBG 

COMMON// G0RHSF(5t20j 65) rHDFPTS(100>100) f ITBLKSZ 
COMMON/BXYCOM/ BXYAR ( 1 0 » 1 0 r 55 > > IJSTGIN(210) »LCND0EC(20) » IBRECr 
1 KFSTGfKPLSTG(70) 

DIMENSION LVECT0R(20) 

C THE THIRD DIMENSION IN THE FOLLOWING THREE ARRAYS MUST EQUAL 
C 'ITBLKSZ'. 

DIMENSION HLAT(5»20,17) »6FLD(5r20»17) »GNBNG(5»20»17) 

EQUIVALENCE < HLAT ( 1 » 1 » 1 ) » 60RHSF ( 1 , 1 , 1 ) ) , ( GFLD ( 1 , 1 » 1 ) f GORHSF ( 1 f 1 » 

1 18) ) f (GNBNG<1» Irl) »G0RHSF(lrl»35) ) 

C 

C USING MDFCDOUBLE FAULT DID CAUSE FAILURE) » THIS FUNCTION 
C COMPUTES THE RATE AT WHICH THE SYSTEM FAILS BY TIME T DUE 
C A CRITICAL FAULT CONDITION, NOTE: THIS FUNCTION IS NOT CALLED 
C UNLESS LVECTOR CONTAINS AT LEAST 2 FAULTS. 

C 

C FAC = 0.0 AT TIME 0. 

FAC =0.0 

C IF 'IBREC' EQUALS 0 OR THIS SUBRUN'S FIRST STAGE DOES NOT EQUAL THE 
C CURRENT COUPLED PAIR SUBSET'S FIRST STAGEr NO CRITICAL PAIR FAULTS 
C EXIST AND THIS FUNCTION EQUALS ZERO. 

RATE WHICH SYSTEMS HAVING L 
FAULTS FAIL AT TIME t DUE TO 
CRITICAL FAULT CONDITIONS 


a' (t|i) 








EXTERNALS 

BUFFIN 

FBCRTL 


TYPE 

ARGS 

4 

REAL 

6 


INLINE FUNCTIONS TYPE 
MOD INTEGER 


ARGS 

2 INTRIN 
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FUNCTION FAPC< ITiLVECTOR) 

C 

COMMON/CONFIG/ NCONUEC < 20 ) » MSRUUEC ( 20 ) , NFLTCAT ( 20 > » MFLTYPE ( 5 » 20 ) » 
1 NSTGS» IFSTG»NTYPS»TRNSFC(5»20) fMAXSTGf MAXCAT,MAXTYP»MXSBRN 
COHMON/CHICOM/ MDP * MEAT » HELD » MBNG » MNBNG » MDF » CHI DBG 
COMMON/NONLDEP/ AXI AR ( 5 f 20 f 65 ) » RXAR ( 20 » 65 ) 

LOGICAL TRNSFC»CHIDBG 

COMMON// G0RHSF(5»20»65) »HDFPTS(100» 100) » ITBLKSZ 
DIMENSION LVECT0R(20) 

C THE THIRD DIMENSION IN THE FOLLOWING THREE ARRAYS MUST EQUAL 
C 'ITBLKSZ'. 

DIMENSION HLAT(5»20r 17) » 6FLD ( 5 » 20 » 1 7 ) » GNBNG ( 5 » 20 » 1 7 ) 

EQUIVALENCE (HLAT< 1 » 1 » 1 ) »GORHSF( 1 » 1 » 1 ) ) » (GFLD( 1 » 1 , 1 ) »GORHSF( 1 » 1 1- 
1 18)) f (GNBNGdrlf 1) fG0RHSF<l»lF35) ) 

THIS FUNCTION COMPUTES THE RATE AT WHICH A SYSTEM HAVING LVECTOR 
FAULTS FAILS AT TIME T DUE TO A NONCRITICAL FAULT. 


RATE AT WHICH SYSTEMS HAVING 
A FAULTS FAIL AT TIME t DUE TO 
ERROR PROPAGATION 

a' (tin = S 

1-r (t) 

PERMANENT 


TRANSIENT 


INLINE FUNCTIONS TYPE ARGS 

mod INTEGER 2 INTRIN 
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FUNCTION FBCRTL( IT» IST6» ICAT» JSTG» JCATfLVECTOR) 

C 

COMMON/CONFIG/ NCONVEC ( 20 ) » MSRVUEC ( 20 ) » NFLTCAT ( 20 > » MFLTYPE ( 5 » 20 ) » 
1 NSTGS»IFSTG»NTYPS>TRNSFC(5»20) »MAXSTG>MAXCAT»MAXTYP»MXSBRN 
COMMON/BXYCOM/ BXYAR ( 1 0 r 10 » 55 ) » IJSTGIN ( 21 0 ) r LCNDVEC ( 20 ) f IBREC > 

1 KFSTG»KPLSTG(70) 

COMMON/NONLDEP/ AXI AR < 5 » 20 » 65 ) » RXAR < 20 » 65 ) 

LOGICAL TRNSFCtTRNSXI »TRNSYJ 
DIMENSION LVECT0R(20) 

DATA MAXFLT/10/ 

C 

C THIS FUNCTION COMPUTES THE PROBABILITY THAT A SYSTEM CONTAINING 
C LVECTOR FAULTS WOULD ENTER AN XlrYJ CRITICAL STATE» REPRESENTED 
C BY ' ISTG» ICAT; JSTGf JCAT' » WERE A CATEGORY YJ FAULT TO OCCUR AT 
C TIME T» REPRESENTED BY 'IT'. 

C 


EXPECTED NUMBER OF 
CRITICAL FAULTS AT TIME t 
GIVEN A PERMANENT FAULTS 




-= / j b (A -VI )p ii„) 

i«i* X/Y X X y x ' x 




P(yy/tUy)C(x^,yj)a^^(t)ay (t) 


Jd S-Zrored //? ( 10^ 10^ sd) 


EXTERNALS 

BXYC 

FPMUX 


TYPE ARGS 


REAL 


•for JeU/)l4iC3n . 
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FUNCTION 


MATHEflATICAL EXPRESSION 



aj^(t)ay(t) 


^i^Yj = PERMANENT 

X 7^ y 


a 2 {t) 


^i/Yj = PERMANENT 

X = y 


IJ (n -£ ) 
X V y' 


ax(t) 


= PERMANENT 
Yj = TRANSIENT 



= TRANSIENT 
Vj = PERMANENT 


^^x" x^ ^ V^y^ = transient 

X 7^ y 

^^x"V^”x"V^^ ^i'^j = transient 

X = y 
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FUNCTION FCLAM( IT» ITAU»ISTG» ICAT) 

C 

COMMON/RATES/ OMGA ( 5 » 20 ) » RLAM ( 5 r 20 > 

COMMON/STEPCOH/ I TSTPS » HAXSTP » RELSTEP > XBASE f NSTGRN f KUT » PSTRNC > 
1 PRCODEfRLPLOTf IAXSRLfSYSMNTfTBCF 

LOGICAL SYSHNTfRLPLOT 
C 


A„ (t) = X t“^i 

Xi 


mr\a 


LirfS c^nr} pc/ie/1'ls E^Qaa-/i'c 


RL^hn(xc/}r^x^^:^) 
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FUNCTION FCYJ( IT, JSTG» JCAT»LVECTOR) 

C 

COMHON/CONFIB/ NCONVEC ( 20 ) , HSRVVEC ( 20 ) , NFLTCAT < 20 ) » MFLTYPE ( 5 , 20 ) r 
1 NSTBSf IFSTG,NTYPS»TRNSFC(5»20) »MAXSTGf MAXCAT,HAXTYP,MXSBRN 
COHMON/CHICOM/ MDP » MLAT , MFLD » MBNG , MNBN6 , MDF , CHI DBG 
LOGICAL TRNSFCrCHIDBG 

COMMON// G0RHBF(5f 20, 65) ,HDFPTS<100,100) , ITBLKSZ 
COMMON/BXYCOM/ BXYAR < 10 , 10 , 55) , IJSTGIN ( 210 ) , LCNDVEC ( 20 ) , IBREC , 

1 KFSTG,KPLSTG(70) 

DIMENSION LVECT0R(20) 

C THE THIRD DIMENTION IN THE FOLLOWING THREE ARRAYS MUST EQUAL 
C 'ITBLKSZ'. 

DIMENSION HLAT(5,20,17) ,GFLD(5,20,17) , GNBNG < 5 , 20 , 1 7 ) 

EQUIVALENCE (HLAT( 1 ,1,1), GORHSF ( 1 , 1 , 1 ) ) , ( GFLD (1,1,1), GORHSF (1,1, 

1 18) ) , (GNBNG ( 1 ,1,1) , GORHSF ( 1, 1 ,35) ) 

C 


PROBABILITY THAT A CATEGORY Yj 
FAULT WOULD PRODUCE A SYSTEM 
FAILURE AT TIME t GIVEN 1 
FAULTS AT TIME f 


C (t|S,) 






(t\i) 


EXTERNALS TYPE ARGS 

FDSCRTL REAL 6 


INLINE FUNCTIONS TYPE 
MOD INTEGER 


ARGS 

2 INTRIN 
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FUNCTION FDSCRTLdTj ISTG,ICAT» JSTG» JCAT»LVECTOR) 

C 

COMHON/CONFIG/ NCONVEC ( 20 ) » MSRVVEC < 20 ) » NFLTCAT ( 20 ) » MFLTYPE ( 5 r 20) » 
1 NSTGSf IFSTG»NTYPS»TRNSFC<5»20) »MAXSTG»MAXCAT»HAXTYP»hXSBRN 
LOGICAL TRNSFC 

COMMON/NONLDEP/ AXI AR ( 5 » 20 » 65 ) » RXAR ( 20 r 65 ) 

COMHON/BXYCOM/ BXYAR (1 0 » 1 0 » 55 ) » IJSTG I N ( 2 1 0 ) » LCNDVEC ( 20 ) dBREC^ 

1 KFSTG»KPLSTG(70) 

DIMENSION LVECT0R(20) 

DATA MAXFLT/10/ 

C 

C THIS FUNCTION COMPUTES THE PROBABILITY THAT A SYSTEM CONTAINING 
C LVECTOR FAULTS WOULD BE IN A SUPERCRITICAL STATE WERE A CATEGORY 
C YJ FAULT» REPRESENTED BY 'JSTG»JCAT'» TO OCCUR AT TIME T, 

C REPRESENTED BY ' IT' . 

C 


I 


EXPECTED NUMBER OF x^y- 
CRITICAL FAULTS^ GIVEN 1 
PERMANENT FAULTS^ THAT WOULD 
BE CREATED AS THE RESULT OF 
A STAGE y FAULT AT TIME t 





(t|i) 



x' y 


X,y 


(A -y 

X 


x' 


A -y )P(y , 
y y 


tiy 


P(Uy,t|jy) 

^x = PERMANENT 

^V^x^^x. = transient 


EXTERNALS TYPE ARGS 

BXYC 5 

FPMUX REAL 4 
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n n fj n o 


FUNCTION FFDFSKMCHI f KINDX»ISTG» ICAT» JSTG» JCAT) 

C 

COMMON/CONFIG/ NCONVEC ( 20 ) » MSRM VEC ( 20 ) , NFLTCAT ( 20 ) » MFLT YPE ( 5 » 20 ) r 
1 NSTGSt IFSTGj NTYPS»TRNSFC( 5»20) >MAXSTGf MAXCAT»MAXTYP»MXSBRN 
COMMON/CHI COM/ MDP » MEAT » MFLD » MBNG » MNBNG » MDF » CH I DBG 
COMMON/NONLDEP/ AX I AR < 5 » 20 t 65 ) » RXAR ( 20 » 65 ) 

LOGICAL TRNSFC»CHIDBG 

COMMON// 60RHSF(5>20»65) »HDFPTS(100»100) j ITBLKSZ 
DIMENSION HBNG<5»20»65) 

EQUIVALENCE (HBNG< 1 » 1 » 1 ) »G0RHSF(1 » 1 » 1 ) ) 

THIS FUNCTION IN COMBINATION WITH SUBROUTINE ABCST REPRESENTS 
SLOW VARYING RELIABILITY FUNCTIONS BY A THIRD-ORDER POLYNOMIAL OVER 
THE RANGE OF INTEREST FOR DOUBLE FAULTS. 


EXTERNALS 

FFSFST 


TYPE ARGS 
REAL 6 
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CJ u o o o 


FUNCTION FFSFSKMCHI >KINDX»ISTG»ICATf JDUMS> JDUMC) 

C 

COMMON/CONFIG/ NCONVEC ( 20 ) » MSRVVEC ( 20 ) » NFLTCAT ( 20 ) , MFLTYPE ( 5 , 20 ) , 
1 NSTGS» IFST6»NTYPS»TRNSFC(5»20) ,MAXSTG»MAXCAT,HAXTYPf MXSBRN 
COMMON/CHICOM/ MDP » MLAT » MFLD » HBNG » MNBNG » HDF » CHI DBG 
COMMON/NONLDEP/ AXI AR < 5 r 20 r 65 ) » RXAR ( 20 y 65 ) 

LOGICAL TRNSFC»CHIDBG 

THIS FUNCTION IN COMBINATION WITH SUBROUTINE ABCST REPRESENTS 
SLOW VARYING RELIABILITY FUNCTIONS BY A THIRD-ORDER POLYNOMIAL OVER 
THE RANGE OF INTEREST FOR SINGLE FAULTS. 


EXTERNALS TYPE ARGS 

FLAM REAL 3 
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FUNCTION FGSKMCHI » IT» ISTG, ICAT) 

C 

COMMON/CHICOH/ MDP > MLAT » MFLD » MBNG » MNBNG » MDF » CHI DBG 
COMMON/CONFIG/ NCONVEC ( 20 ) » MSRUVEC ( 20 ) » NFLTCAT ( 20 ) » HFLTYPE ( 5 » 20 ) » 
1 NSTGS> IFSTG»NTYPSjTRNSFC< 5»20) rMAXST6»MAXCAT»MAXTYP»MXSBRN 
LOGICAL TRNSFC jCHIDBG 

COHMON/NONLDEP/ AXI AR ( 5 » 20 » 65 ) » RXAR (20 » 65 ) 

COMMON// G0RHSF(5»20»65) » HDFPTS( 100 , 100 > f ITBLKSZ 
C 

' C THIS FUNCTION COMPUTES THE PROBABILITY THAT A CATEGORY XI FAULT IS 
C IN STATE MCHI AT TIME T» REPRESENTED BY IT> GIVEN THAT IT WAS LATENT 
C PREVIOUSLY. 
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FUNCTION FHDFSKMCHI . I T , I STG » ICAT » JSTG » JC AT ) 

C 

COMMON/STEPCOM/ I TSTPS » MAXSTP » RELSTEP » TBASE » NSTGRN > KUT » PSTRNC , 

1 PRCODE»RLPLOT»IAXSRL»SYSMNT»lBCF 

COHHON/CONFIG/ NCONVEC < 20 ) » MSRVVEC ( 20 ) » NFLTCAT ( 20 ) » MFLTYPE ( 5 » 20 ) » 

1 NSTGS»IFSTG»NTYPS»TRNSFC(5»20) » MAXSTG » HAXCAT » MAXT YP » MXSBRN 
COMMON/CHICOM/ HOP i MLAT > HFLD » MBNG » HNBNG » MDF , CH I DBG 
COMMON/CVRGCOM/ CHSTO < 5 » 65 » 5 > » CMST 1 ( 5 » 65 » 5 > » CMST2 < 5 » 65 > 5 ) » 

1 CHDF0(5»5»65> »CMDF1 (5»5»65) » CMHr2 ( 5 » 5 » 65 ) f 

2 TZER0ST(5»5) »T0DF(5»5> »TRANS(5> 

LOGICAL TRANS fTRNSFC»CHIDB6»SYSMNT»RLPL0T 
EXTERNAL FFDFST 

C 

C THIS ROUTINE COMBINES A POLYNOMIAL REPRESENTATION OF RELIABILITY 
C FUNCTIONS(SLOU VARYING) AND MOMENTS OF COVERAGE PROBABILITY FUNCTIONS 
C (FAST VARYING) TO YIELD A CONVOLUTED INTEGRATION OF THE TWO TYPES OF 
C DOUBLE FAULT FUNCTIONS. 

C 


RATE AT WHICH AN -CRITICAL 
FAULT CAUSES SYSTEM FAILURE 

hpp(tlxi,yj) = apj,(t|Xj^,yj)m^p(t|x.,yj) 

+ bpp(tlx.,y.)mJp(t|xj^,yj) 

+ Cpp(t|Xi,y.)m^p(t|Xi,yj) 


EXTERNALS 

TYPE 

ARGS 

ABCST 


11 

FFDFST 


0 

INLINE FUNCTIONS 

TYPE 

ARGS 

MOD 

INTEGER 

2 INTRIN 
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noon non 


FUNCTION FHSFSKMCHI » I T , I ST6 » ICAT ) 

C 

COHMON/STEPCOM/ I TSTPS f MAXSTP » RELSTEP » XBASE » NSTGRN » KWT » PSTRNC » 

1 PRCODE>RLPLOT»IAXSRL»SYSMNTf TBCF 

COMHON/CONFIG/ NCONVEC ( 20 ) , MSRVUEC ( 20 ) » NFLTCAT ( 20 ) , HFLTYPE ( 5 , 20 ) , 

1 NSTGS» IFSTG»NTYPS»TRNSFC<5»20> »HAXSTG»MAXCATf MAXTYP»MXSBRN 
COHMON/CHICOM/ MDP » MLAT f MFLD » MBNG f MNBNG » MDF f CHI DBG 
COMMON/CVRGCOH/ CMSTO ( 5 » 65 » 5 ) » CMSTl ( 5 » 65 » 5 ) , CMST2 ( 5 » 65 » 5 ) f 

1 CMDF0<5»5.65),CMDF1 (5»5»65> f CHDF2(5»5»65> T 

2 TZEROST<5»5)»TODF(5»5) fTRANS(5) 

LOGICAL TRANS »TRNSFC» CHI DBG »SYSMNT»RLPLOT 
EXTERNAL FFSFST 
DATA JDUM/0/ 

THIS ROUTINE COMBINES A POLYNOMIAL REPRESENTATION OF RELIABILITY 
FUNCTIONS(SLOW VARYING) AND MOMENTS OF COVERAGE PROBABILITY FUNCTIONS 
(FAST VARYING) TO YIELD A CONVOLUTED INTEGRATION OF THE TWO TYPES OF 
SINGLE FAULT FUNCTIONS. 


EXTERNALS 

ABCST 

FFSFST 


TYPE 


ARGS 

11 

0 


INLINE FUNCTIONS TYPE 
MOD INTEGER 


-ARGS 

2 INTRIN 


See the next page for the Single-fault H functions computed using 
this function. 
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DEFINITION 


FUNCTION MATHEMATICAL EXPRESSION 


PROBABILITY THAT A CATEGORY x^- 
TRANSIENT FAULT IS DETECTED AS 

PERMANENT 

Hjjp^(tlx^) = aj^p(t|x^)m^p(tlx^)+bj^p(t|x^)mjp(tlx^) 

+ Cpp{t|x.)xn 2 p(t|x.) 

PROBABILITY OF A LATENT CATEGORY 
FAULT AT TIME t 


HL(tlx.) 

RATE OF ERROR PROPAGATION 
FAILURE DUE TO A CATEGORY x^ 
FAULT 

hp(t|x^) = 


PROBABILITY OF A BENIGN 
LATENT FAULT AT TIME t 

HgCtlx.) 


PROBABILITY OF A NON-BENIGN 
LATENT FAULT AT TIME t 

Hpitlx.) 


aL(tlXi)M°(tlXi)+bL(tlXi)Mj;(tlXi) 

ap(t|x^)m^{t|x^)+bp(t|x^)mp(tlx^) 
+ Cp(t]x^)m^(tlx^) 

a3(tlXi)M°(tlXi)+bB(tlXi)Mg(t|Xj.) 
+ CB(t|x^)MB(t|x^) 

a^(t|Xi)M^(tIxi)+b-(tlXi)M^(t|xj^) 
+ Cgltlx-lMl (t|x.) 
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FUNCTION FINT6RT( IT. ARTOINT) 

C 

COMMON/STEPCOM/ I TSTPS » MAXSTP » RELSTEP » TBASE » NSTGRN » KWT » PSTRNC . 
1 PRCODE.RLPLOT, lAXSRL.SYSMNT.TBCF 

DIMENSION ART0INT(3) 

LOGICAL SYSMNT.RLPLOT 
C 

C THIS FUNCTION USES THE 1/3 SIMPSON'S INTEGRATION METHOD ON THE 
C VALUES STORED IN ARTOINT WHEN *IT’.GT.2. 'IT' REPRESENTS 
C WHICH TIME STEP IS BEING COMPUTED. 

C WHEN 'IT* EQUALS 2. THE TRAPEZOIDAL INTEGRATION METHOD IS USED. 

C 
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FUNCTION FLAMdTrlSTGrlCAT) 

C 

COMMON/RATES/ 0MGA(5»2O) »RLAM(5»20) 

COMMON/STEPCOH/ I TSTPS t HAXSTP » RELSTEP » TBASE » NSTGRN i KWT » PSTRNC f 
1 PRCODEfRLPLOT i I AXSRL f SYSHNT » TBCF 

LOGICAL SYSMNTrRLPLOT 
C 


X (t) 


<■> 


X. 



t 



L)anal?le? 




RLfhV xStz^) 

0 ) 016 /} ^c.fh^ ^ 
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noon 


FUNCTION FNFACT(N) 

THIS FUNCTION COMPUTES N FACTORIAL 
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o o n n n 


FUNCTION FPMUX(IT»IST6iHUX»LX) 

C 

COMMON/CONFIG/ NCONVEC ( 20 ) » MSRVVEC ( 20 ) , NFLTCAT ( 20 ) » MFLTYPE ( 5 » 20 ) ^ 
1 NSTGS» IFST6»NTYPS,TRNSFC(5r20) »MAXSTG»MAXCAT»MAXTYP»MXSBRN 
COHMON/NONLDEFV AXI AR ( 5 » 20 » 65 ) » RX AR ( 20 , 65 ) 

LOGICAL TRNSFC 

THIS FUNCTION COMPUTES THE PROBABILITY THAT A SUBSYSTEM CONTAINS 
MUX STAGE X LATENT FAULTS GIVEN THAT IT HAS LX FAULTY STAGE X 
MODULES. 


PROBABILITY THAT A SUBSYSTEM 
CONTAINS STAGE ^ LATENT 
PERMANENT FAULTS GIVEN THAT 
IT HAS STAGE x PERMANENT 
FAULTS 







X X 



EXTERNALS 

FNCK 


TYPE ARGS 
REAL 2 
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n n o n 


FUNCTION FPSTAR( ITfLVECTOR) 

C 

COMMON/CONFIG/ NC0NVEC(20) » MSRUVEC ( 20 ) , NFLTCAl ( 20 ) » MFLT YPE ( 5 » 20 ) » 
1 NSTGSf IFST6>NTYPS»TRNSFC(5»20) » H AXSTG » MAXCAT » HAXT YP > MXSBRN 
COHMQN/NDNLDEP/ AXI AR ( 5 » 20 » 65 ) > RXAR ( 20 » 65 ) 

LOGICAL TRNSFC 

THIS FUNCTION COMPUTES THE PROBABILITY THAT A SYSTEM HAS SUSTAINED 
EXACTLY LVECTOR FAILURES BY TIME T» REPRESENTED BY IT. 

DIMENSION LVECT0R(20) 


EXTERNALS 

FNCK 


type ARGS 
REAL 2 
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n n n n n 


FUNCTION FPSTREC< IT» JSTG»LVECTOR»PSTPRV> 

C 

COHMON/CONFIG/ NCONVEC ( 20 ) » MSRVUEC ( 20 ) » NFLTCAT ( 20 ) r HFLT YPE ( 5 » 20 ) » 
1 NSTGS» IFSTG»NTYPS»TRNSFC(5,20) ,MAXSTG»MAXCAT»MAXTYP»MXSBRN 
COMMON/NONLDEP/ AXI AR < 5 » 20 » 65 ) » RXAR ( 20 » 65 ) 

LOGICAL TRNSFC 

THIS FUNCTION COMPUTES PSTAR FOR LVECTOR RECURSIVELY - GIVEN 
PSTAR FOR LVECTOR MINUS THE UNIT VECTOR REPRESENTED BY JSTG 
PASSED AS PSTPRV. 

DIMENSION LVECT0R(20) 
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n n n o n 


FUNCTION FRXIFF( IT» ISTGfICATfTRNSFLT) 

C 

COMMON/CONFIG/ NCONVEC ( 20 ) » MSRVVEC ( 20 ) i NFLTCAT ( 20 ) » MFLTYPE ( 5 f 20 ) , 

1 NSTGSj IFSTG»NTYPS»TRNSFC( 5f20) » MAXSTG » HAXCAT r MAXT YP » MXSBRN 
COMMON/CHICOM/ MDP » MEAT » MFLD , MBNG r MNBNG » MDF , CHI HBG 
COMMON/CVRGCGH/ CHSTO ( 5 » 65 » 5 ) » CMST 1 ( 5 » 65 » 5 ) » CMST2 ( 5 » 65 » 5 ) , 

1 CHDF0(5»5»65) »CMDF1 (5»5»65) » CMDF2 ( 5 f 5 » 65 ) » 

2 TZER0ST<5»5)»T0DF(5»5) »TRANS(5) 

COMMON/NONLDEP/ AXI AR < 5 » 20 » 65 ) f RXAR ( 20 » 65 ) 

COMMON// G0RHSF<5f20f 65) fHDFPTSdOOf 100) » ITBLKSZ 
LOGICAL TRANS rTRNSFC» CHI DBG »TRNSFLT 
DIMENSION HDPFC(5»20»3) fHINT(5»20r3) f HDP(3) 

C USE AXIAR AS WORKING STORAGE FOR THIS ROUTINE BECAUSE AXIAR IS NOT 
C COMPUTED UNTIL RXAR IS COMPLETELY DEFINED, 

EQUIVALENCE ( AXIAR ( 1 d » 1 )» HDPFC ( 1 » 1 > 1 ))»( AXIAR ( 1 r 1 , 4 ), HINT ( 1 » 1 d ) ) 

THIS FUNCTION COMPUTES THE PROBABILITY THAT A GIVEN STAGE X 
MODULE HAS EITHER NOT EXPERIENCED» OR HAS RECOVERED FROM» ANY 
CATEGORY XI FAULTS BY TIME T. 


PROBABILITY THAT A GIVEN STAGE 
X MODULE HAS NOT EXPERIENCED A 

category-xjl fault by time t 


I ~^x- 

I ® ^ PERMANENT 

i 

TRANSIENT 


EXTERNALS TYPE ARGS 


FCLAM REAL 4 

FHSFST REAL 4 

FINTGRT REAL _2 

PREEXP REAL 1 
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3 .4 General Purpose Routines 

The following routines are used in at least two of the three 
program modules described previously. NOTE: Subroutine CAREPLT 

interfaces to NASA's Graphic Plotting Package, and subroutine 
CPLOT interfaces to the DISSPLA Plotting Package. Currently, 
programs CVGPLT and RELPLT call CPLOT. Slight modifications must 
be made to these programs in order to use CAREPLT. 
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3.4.1 General Purpose Siibroutines 


SUBROUTINE BUFBLK ( I NOUT » lUNIT » BLOCK » ILST » FATAL » E0FFL6 > 

C 

C THIS SUBROUTINE BUFFERS IN A BLOCK OF DATA INTO ' BLOCK ' (1 >.♦' BLOCK ' 
C (ILST)» IF 'INOUT' ,EQ.1» FROM UNIT 'lUNIT'. IT BUFFERS OUT A BLOCK 
C OF DATA FROM ' BLOCK '( 1 BLOCK '( ILST ) r IF ' I NOUT ' . EQ . 2 » TO UNIT 
C 'lUNIT'. 

C 

DIMENSION BLOCK(l) 

LOGICAL FATALrEOFFLG 


EXTERNALS 

UNIT 


TYPE ARGS 
REAL 1 
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SUBROUTINE CAREPLT ( ARPASD » XPASD » STEP , 6ENXPTS r NPTS » LN0RL6 ) 

C 1l.1ti**********tt*t**1lit***ttt*tt**1lit***1li**t*t**t*t*t*t**1li*illiilli1li**1li*t*t*tt* 
C THIS ROUTINE WORKS WITH NASA'S GRAPHIC PLOTTING PACKAGE, 

^b ^b ^b ^b ^b ^b ^b ^ ^b ^b ^b 

^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ * 

DIMENSION ART0PLT<515) »X(515) »ARPASD(1) fXPASDd) 

C*****NOTE THAT TWO ADDITIONAL ELEMENTS ARE ASSIGNED AS 
C****«REQUIRED BY ASCALE » LINPLT > SCALOG 
LOGICAL GENXPTS 

DATA XAXlSMX/8.0/» YAXISMX/10.0/»XDV/2.0/f XTIC/.5/» 

1 YDV/1 ,0/»YTIC/l .0/ 

C 

^|c ^|c 4^ 4^ 4^ 4^ 4^ 4^ 4^ 4^ 4^ ^|c 4^ 4^ 4^ 4^ 2^ 2|( 2^ 4^ 2^ 2^ 4^ 2§C 2^ 2^ 2|C 2^ 2^ 

c 

c 

C ARTOPLT 

C THIS IS THE ARRAY TO BE PLOTTED WITH ’STEP* AS STEP SIZE 

C NPTS 

C THIS IS THE NUMBER OF POINTS TO BE PLOTTED 

C CAN BE PLOTTED, 

C lAXSTYP 

C THIS IS THE Y-AXIS TYPE DESIRED ( LINEAR » SEMI-LOG OR BOTH) 

C XAXISMX 

C THIS IS THE LENGTH IN INCHES OVER WHICH THE RANGE 

C OF X VALUES IS TO BE PLOTTED. 

C XDV 

C THE NUMBER OF DIVISIONS PER INCH OF PAPER TO BE 

C USED FOR THE X AXIS. 

C XTIC 

C THIS IS THE LENTGTH IN INCHES BETWEEN ANNOTATED 

C TIC MARKS ON THE X AXIS. 

C YAXISMX 

C THIS IS THE LENGTH IN INCHES OVER WHICH THE RANGE 

C OF Y VALUES IS TO BE PLOTTED. 

C YDV 

C THE NUMBER OF DIVISIONS PER INCH OF PAPER TO BE 

C USED FOR THE Y AXIS. 

C YTIC 

C THIS IS THE LENGTH IN INCHES BETWEEN ANNOTATED 

C TIC MARKS ON THE Y AXIS 


LS 

TYPE 

ARGS 

ASCALE 


5 

AXES 


11 

AXESLOG 


10 

CALPLT 


3 

LGLIN 


8 

LINPLT 


8 

NFRAME 


0 

PSEUDO 


0 

SCALOG 


5 
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SUBROUTINE CPLOT ( ARTOPLT i X » STEP f DENXPTS » NPTS » LNORLG » IT ITLE » 

1 JTITLE) 

C aia llf ^ir 

^p ^p ^k ^p ^p ^p 4^ 4^ 4^ ^k 4^ ^k ^k ^k ^p 

C THIS ROUTINE WORKS WITH THE DISSPLA PLOTTING PACKAGE, 

C ^i ^lf ^ ^k ^k ^k ^k ^k ^k ^k ^k ^k ^k ^k ^k ^k ^k ^k ^k ^k ^k 

^p ^p ^p ^p ^p ^p ^p ^p ^p ^p 4^ ^p 4^ 4^ ^p ^p ^p ^P ^P ^k ^p ^P ^p ^P ^P ^P ^P ^P ^P ^P ^P ^P ^P 

C NOTE - ARRAYS 'ARTOPLT' AND 'X' MUST BE OF THE SAME DIMENSION 
C IN THE CALLING ROUTINE, 

DIMENSION ARTOPLT(l) »X(1 ) fITITLE(I) r JTITLE(l) 

LOGICAL GENXPTS 


EXTERNALS 

ALGPLT 

AXSPLT 

BGNPL 

BLNKl 

COMPRS 

CURVE 

ENDPL 

FRAME 

GRAPH 

GRID 

MESSAG 

NOCHEK 

RESET 

TITLE 

YLOG 


TYPE ARGS 

5 

6 
1 
5 
0 
A 
1 
0 
4 

4 

0 

1 

8 

4 
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o n n n n 


3.4.2 General Purpose Functions 


FUNCTION FNCK<NFAC»KFAC) 

THIS FUNCTION COMPUTES BINOMIAL COEFFICIENTS. 

GIVEN TWO PARAMETERS N»K WILL COMPUTE NCR I.E. NFACT/ ( N-K ) FACT 
* KFACT. THIS IS EQUIVALENT TO N ( N- 1 )...( K+ 1 >/( N-K ) FACT . 
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FUNCTION PREEXP<X) 


DATA REALMAX/1 .0E+322/»REALMIN/l .0E-293/»EXPMAX/741 ,67/f 
1 EXPHIN/-675.82/ 


IF(X.GT.EXPMIN .AND* X.LT.EXPMAX) GO TO 100 
C SET FUNCTION TO A VALUE VERY CLOSE TO 0.0 BUT NOT EQUAL TO 0.0 
IF(X.LE.EXPMIN) PREEXP = REALMIN 
C SET FUNCTION TO THE MAXIMUM NUMBER THE CDC CAN HANDLE. 
IF(X.GE.EXPMAX) PREEXP = REALMAX 
GO TO 200 

100 PREEXP = EXP<X) 

200 RETURN 
END 


EXTERNALS 

EXP 


TYPE ARGS 
REAL 1 LIBRARY 
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4 .0 Error Messages 


The following sections contain the error and warning messages 
generated by modules CAREIN, COVRGE, and CARES. The utility rou- 
tine error messages can occur in any of the three modules . The 
error messages are fatal to the computer run; therefore, the run 
will halt after one occurs. The warning messages are nonfatal and 
are used to alert the user to any out-of-the-ordinary processing. 

The variables enclosed in single quote marks are used to show 
which variables will have their contents printed. The following 
variables, used in the error messages, have values which do not 
change during the course of the computer run: 


Variable Name 

Value 

Description 

MAXTYP 

5 

Maximum Niamber of 
Fault Types 

MXSTGS 

70 

Maximum Number of 
Stages 

MAXCAT 

5 

Maximum Number of 
Fault Categories 

MXUNTS 

70 

Maximum Number of 
Coupled Units Per 
Fault Tree 

MXCPLS 

20 

Maximum Number of 
Coupled Stages Per 
Fault Tree 

MXSBRN 

35 

Maximum Number of 
Subruns 

INMAX 

513 

Maximum Index Into 
Coverage Functions 

MXSTGF 

13 

Maximum Number of 
Stage Failures 
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Variable Name 

Value 

Description 

UMAX 

210 

Maximum Number of 
Unique Fault Pair 
Functions Per 
BXYFL Subfile 

MAXSTG 

20 

Maximum Number of 
Stages Per Subrun 


Beside each message is a code used to differentiate the user 
errors from the possible internal program errors for which tests 
are made . 
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4.1 CAREIN ERROR MESSAGES (FATAL) 


** ERROR - NFTYPS = 'NFTYPS' MUST BE LESS THAN OR EQUAL TO 'MAXTYP' AND UE 

NSTGES = 'NSTGES' MUST BE LESS THAN OR EQUAL TO 'MXSTGS'. 

** ERROR - NFCATS ( ' ISTG ' ) = 'NFCATS (ISTG) ' MUST BE LESS THAN OR EQUAL TO 'MAXCAT'. UE 

** ERROR - EOF ENCOUNTERED ON FILE CREIN WHEN MORE DATA WAS REQUIRED. UE 

** ERROR - STAGE NUMBERS MUST RANGE FROM 1 TO 'MXSTGS' MAXIMUM UE 

AND THE LAST STAGE MUST BE LESS THAN OR EQUAL TO 'NSTGES' DEFINED AS 'NSTGES' 

- NOT 'IFSTG' TO ' ILSTG ' . 

** ERROR - EXCEEDED MAXIMUM NUMBER OF COUPLED UNITS UE 

PER CRITICAL PAIR FAULT-TREE OF 'MXUNTS'. 

** ERROR - ILLEGAL STAGE SPECIFIED IN CRITICAL PAIR DATA = 'ISTG'. UE 

** ERROR - NUMBER OF UNITS SPECIFIED FOR STAGE 'ISTG', IN THE CRITICAL PAIR DATA, UE 

DOES NOT EQUAL THE NUMBER OF IDENTICAL L’NITS SPECIFIED FOR THIS STAGE. 

IF ALL UNITS ARE NOT INVOLVED IN THE CRITICAL PAIR FAULT-TREE, SET THE NUMBER 
OF OPERATING STATES TO THE NUMBER OF INVOLVED UNITS FOR THIS STAGE, AND 
INCLUDE EACH UNIT IN THE FAULT-TREE . 

** ERROR - EXCEEDED MAXIMUM NUMBER OF COUPLED STAGES UE 

PER CRITICAL PAIR FAULT-TREE OF 'MXCPLS'. 

** ERROR - EXCEEDED MAXIMUM RANGE OF COUPLED STAGES UE 

PER CRITICAL PAIR FAULT-TREE OF 'MXCPLS', 

IN THE TREE HAVING FIRST STAGE = 'KSTGF' AND THE LAST STAGE = 'KSTGL' 


NOTE: UE = USER ERROR; IE = INTERNAL ERROR. 
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CARE IN ERROR MESSAGES (FATAL) 


** ERROR - COUPLED STAGES, DEFINED IN CRITICAL PAIR FAULT -TREE S , UE 

MUST BE IN ASCENDING ORDER (FROM ONE TREE TO THE NEXT ALSO) . 

** ERROR - OCCURRED WHILE PROCESSING MINTERM FILE FOR CRITICAL PAIR FAULT DEFINITIONS. IE 

■PRBMT' ' (MNTRMV(I) ,I=1,NUNTS) ' 

** ERROR - SUBRUN ARRAY INDEX GREATER THAN MAXIMUM NUMBER OF SUBRUNS = 'MXSBRN'. IE 

** ERROR - ILLEGAL FIRST STAGE IN CRITICAL-PAIR SUBSET ARRAY: KSSETS ( 1 , ' ISBSET ' ) . IE 

CAREIN WARNING MESSAGES (NON-FATAL) 

** WARNING - INPUT FILE 'CREIN' EMPTY - DEFAULTS WILL BE USED. 

** WARNING - NUMBER OF INTEGRATION STEPS, PLUS 1 TO INCLUDE TIME ZERO, 

REDUCED FROM 'NSTEPS' TO THE MAXIMUM OF 'MAXSTP ' . 


NOTE : UE 


USER ERROR; IE 


INTERNAL ERROR 
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4.2 COVRGE ERROR MESSAGES (FATAL) 


** ERROR - NUMBER OF FAULT TYPES REQUESTED = 'NTYPS' IE 

BUT MUST BE LESS THAN OR EQUAL TO 'MAXTYP*. 

** ERROR - INCORRECT RELIABILITY TIME BASE = 'TBASE'. IE 

** ERROR - THRASHING OCCURRED IN SUBROUTINE COMPFUN. 

RERUN PROGRAM WITH A DIFFERENT ' DBLDF ' VALUE. 

** ERROR - ITH G SINGLE-FAULT FUNCTION = 'ITH' WHILE ONLY G FUNCTIONS 1 THROUGH 9 IE 

MAY BE COMPUTED WITH FUNCTION FGSNGL. 

** ERROR - 'DELTA* ( 'ITYP') EQUALS ZERO, WHILE THE CONSTANT DENSITY FUNCTION WAS UE 

CHOSEN FOR 'DELTA'. 

EITHER ASSIGN 'DELTA' A VALUE OR SPECIFY CONSTANT RATE FUNCTION FOR THIS FAULT TYPE. 

** ERROR - 'RHO ' ( ' ITYP ' ) EQUALS ZERO, WHILE THE CONSTANT DENSITY FUNCTION WAS CHOSEN UE 
FOR 'RHO'. 

EITHER ASSIGN 'RHO' A VALUE OR SPECIFY CONSTANT RATE FUNCTION FOR THIS FAULT TYPE. 

** ERROR - BOTH 'RHO' AND 'DELTA' ('ITYP') EQUAL ZFRO, WHILE THE CONSTANT DENSITY UE 

FUNCTION WAS CHOSEN FOR BOTH FAULT RATES. 

EITHER ASSIGN THEM A VALUE OR SPECIFY CONSTANT RATE FUNCTIONS FOR THIS FAULT TYPE, 

** ERROR - ITH C DOUBLE-FAULT FUNCTION = 'ITH' IffllLE ONLY C FUNCTIONS 1 AND 2 MAY BE IE 

COMPUTED WITH FUNCTION FCDBL. 

** ERROR - ITH F DOUBLE-FAULT FUNCTION = 'ITH' WHILE ONLY F FUNCTIONS 1 AND 2 MAY BE IE 
COMPUTED WITH FUNCTION FFDBL. 


NOTE : UE 


USER ERROR; IE 


INTERNAL ERROR 
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COVRGE ERROR MESSAGES (FATAL) 


** ERROR - ITH B DOUBLE-FAULT FUNCTION = ’ ITH ' WHILE ONLY B FUNCTIONS 1 AND 2 MAY BE IE 
COMPUTED WITH FUNCTION FBDBL. 

** ERROR - THRASHING OCCURRED IN SUBROUTINE SUMARS . RERUN PROGRAM WITH A DIFFERENT UE 

■DBLDF' VALUE. 

** ERROR - THRASHING OCCURRED IN SUBROUTINE VOLTERA. RERUN PROGRAM WITH A DIFFERENT UE 
'DBLDF' VALUE. 

** ERROR - THRASHING OCCURRED IN SUBROUTINE VLTNREC . RERUN PROGRAM WITH A DIFFERENT UE 
' DBLDF ' VALUE . 

** ERROR - INDEX PASSED TO FUNCTION FTCHSTP = ' INDX ' IS LARGER THAN THE TOTAL NUMBER IE 
OF POINTS COMPUTED. 

** ERROR - ILLEGAL MCHI OF 'MCHI ' WHICH SPECIFIES WHICH SINGLE OR DOUBLE-FAULT IE 

FUNCTION TO STORE INTO COMMON /CVRGCOM/ , 

MUST BE IN THE RANGE OF 1 THROUGH 6 ONLY. 

** ERROR WITH MOMNT = 'MOMNT' IN FUNCTION SIMPINT. IE 

ONLY MOMENTS 0,1, AND 2 ARE VALID. 

** ERROR WITH INDICES REPRESENTING INTEGRATION LIMITS IN FUNCTION SIMPINT: 'ITFROM' = IE 

'ITFROM' 'ITTO' = ' ITTO ' WHERE 'INMAX' = 'INMAX'. 

** ERROR - INDEX INTO 'TMAR' = 'IT' IS GREATER THAN MAXIMUM ALLOWABLE INDEX OF 'INMAX'. IE 


NOTE: UE = USER ERROR; IE = INTERNAL ERROR. 
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COVRGE WARNING MESSAGES (NON-FATAL) 


** WARNING - NSTPIN GREATER THAN 'NSTPMX • IN SUBROUTINE COMPFUN . ZERODF DECREASED 
TO ' ZERODF ' . 

** WARNING - NUMBER OF POINTS REQUIRED TO DEFINE A FUNCTION ARRAY IS LARGER THAN 
' INMAX ' . ZERODF INCREASED TO ' ZERODF ' . 

** WARNING - NSTPIN GREATER THAN 'NSTPMX' IN SUBROUTINE S UMARS . ZERODF DECREASED 
TO ' ZERODF ' . 

** WARNING - NPSTPIN GREATER THAN 'NSTPMX'. ZERODF DECREASED TO 'ZERODF'. 

** WARNING - NUMBER OF POINTS REQUIRED TO DEFINE P ARRAY IS LARGER THAN 'INMAX' . 
ZERODF INCREASED TO 'ZERODF'. 
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4.3 CARE 3 ERROR MESSAGES (FATAL) 


** ERROR - INCORRECT TIME BASE = ' TBASE ' . IE 

** ERROR - 'KWT' = 'KWT'. ' KWT ' MUST BE INPUT WHEN THE SYSTEM MINTERM FILE 'FT15F' UE 

IS MISSING. 

THE NUMBER OF STAGE FAILURES CANNOT EXCEED 'MXSTGF'. 

** ERROR - SUBRUN NUMBER EXCEEDED MAXIMUM OF 'MXSBRN'. IE 

** ERROR - EOF ENCOUNTERED ON UNIT 4 DURING BUFFER IN. IE 

** ERROR - WITH SYSTEM MINTERM FILE 'FT15F' - IE 

'PRBMT' (MNTRMV(ISTG) , ISTG=1 ,NSTGRN) 

** ERROR - NUMBER OF UNIQUE FAULT PAIR FUNCTIONS = 'NBXY' FOR 'NSTGS' COUPLED STAGES, IE 
WHICH EXCEEDS THE MAXIMUM OF 'UMAX' FOR THE MAXIMUM NUMBER OF COUPLED STAGES = 

'MAXSTG' . 

** ERROR - INVALID BXYAR INDEX OF ' INBXY ' EXISTS IN INDEX ARRAY IJSTGIN ( ' I JSTG ' ) . IE 

** ERROR - INVALID CONDITIONAL FAULT VALUE OF ' LC ' EXISTS IN CONDITIONAL FAULT UE 

VECTOR LCNDVEC ( ' I ' ) . 

** ERROR - INVALID SPECIFICATION OF WHICH BXY FUNCTION DEFINITION RECORD IS IE 

CURRENTLY IN MEMORY - IBP^C = ' IBREC ' . 

** ERROR - EOF ENCOUNTERED ON PREVIOUS BUFFER IN ON UNIT 'lUNIT' WHILE TRYING TO IE 

READ RECORD ' IT ' . 


NOTE: UE = USER ERROR; IE = INTERNAL ERROR. 
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CARES ERROR MESSAGES (FATAL) 


** ERROR - PARITY ERROR ENCOUNTERED DURING BUFFER IN ON UNIT 'lUNIT' WHILE TRYING TO IE 
READ RECORD ' IT ' . 

** ERROR - EOF ENCOUNTERED ON PREVIOUS BUFFER OUT ON UNIT 'lUNIT' WHILE TRYING TO IE 

WRITE RECORD ' ITB ' . 

** ERROR - PARITY ERROR ENCOUNTERED DURING BUFFER OUT ON UNIT 'lUNIT' WHILE TRYING IE 

TO WRITE RECORD 'ITB*. 

** ERROR - INVALID SPECIFICATION OF WHICH BXY FUNCTION DEFINITION RECORD IS CURRENTLY IE 
IN MEMORY - IBREC = ' IBREC ' . 

** ERROR - CURRENT BLOCK OF CRITICAL PAIR DATA IN MEMORY, STARTING AT STAGE 'KFSTG' IE 
DOES NOT MATCH THE CURRENT SUBRUN STARTING AT STAGE ' IFSTG ' . 

** ERROR IN FPMUX -'MUX' LATENT FAULTS IN STAGE 'ISTG' IS LARGER THAN THE TOTAL IE 

NUMBER OF FAULTS = 'LX' . 

** ERROR - FGST CALLED WITH AN INVALID MCHI OF 'MCHI'. IE 

** ERROR - INCORRECT FAILED STATE REQUESTED FOR USE IN FUNCTION FHSFST: IE 

MCHI = 'MCHI'. 

** ERROR - ILLEGAL FAULT TYPE = ' ITYP ' WHILE TRYING TO READ COVERAGE FUNCTIONS IN IE 

ROUTINE FHSFST. 

** ERROR - INCORRECT FAILED STATE REQUESTED FOR USE IN FUNCTION FHDFST: IE 

MCHI = 'MCHI'. 


NOTE: UE = USER ERROR; IE = INTERNAL ERROR. 
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CARE 3 ERROR MESSAGES (FATAL) 


** ERROR - ILLEGAL FAULT TYPE (S) = ' ITYP ' AND/OR 'JTYP' WHILE TRYING TO READ IE 

COVERAGE FUNCTIONS IN ROUTINE FHDFST. 

** ERROR - INCORRECT FAILED STATE REQUESTED FOR USE IN FUNCTION FFSFST: IE 

MCHI = 'MCHI'. 

** ERROR - INCORRECT FAILED STATE REQUESTED FOR USE IN FUNCTION FFDFST: IE 

MCHI = 'MCHI'. 

** ERROR - ILLEGAL NUMBER PASSED TO N FACTORIAL FUNCTION = 'N'. IE 


CARE3 WARNING MESSAGES (NON-FATAL) 


** WARNING - SYSTEM MINTERM FILE 'FT15F’, GENERATED BY PROGRAM CAREIN, IS EMPTY. 
CARE3 WILL GENERATE P*'S WITHOUT USING THE SYSTEM FAULT-TREE DEFINITION, 

SINCE A 'KWT' VALUE OF 'KWT' WAS SPECIFIED. 

** WARNING - CRITICAL FAULT PAIR FILE 'BXYFL' IS EMPTY. 

THIS RUN ASSUMES THAT NO CRITICAL PAIRS OF FAULTS EXIST. 


NOTE: UE = USER ERROR; IE = INTERNAL ERROR 
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4.4 UTILITY ROUTINE ERROR MESSAGES (FATAL) 


** ERROR - ILLEGAL COMBINATORIAL; NFAC = 'NFAC KFAC = 'KFAC. IE 
** ERROR - ILLEGAL VALUE FOR ' INOUT ' IN SUBROUTINE BUFBLK . IE 
** ERROR - EOF ENCOUNTERED ON UNIT 'lUNIT' DURING BUFFER 'BTYPE'. IE 
** ERROR - PARITY ERROR ENCOUNTERED ON UNIT 'lUNIT' DURING BUFFER 'BTYPE'. IE 


NOTE: UE = USER ERROR; IE = INTERNAL ERROR. 
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The \/o1terra integral equation has received much attention in 

the numerical analysis literature. Many al 90 ''ithms exist for solvina 

the Volterra equation due to the fact that "there is a st'ong link 

between Volterra equations ano initial value orainary aifferential 

equation p r oo 1 em s . . . ( bu t ) the major difference in procedure/ compared 

with differential equations/ lies in the error Cand step size control) 

tecnniques needed ...( wh i 1 e ) no algorithms exist for "stiff" volterra 
1 

problems." 

The purpose of this paper is not to descrioe existing algorithms 
for Volterra integral equations but rather to explore one type of 
Volterra equation which plays a significant role in reliability theory - 
the renewal equation. The renewal equation is a linear Volterra 
equation with convolution kernal/ defined as the expected number of 
renewals by time t. In particular, the "stiff" renewal equation 
will be explored. This is an equation where the functions involved 
contain vastly oifferent time constants. The computational complexi- 
ties involvea in defining an algorithm to solve this equation will 
be explored; and a solution will be proposed based on the generalized 
trapezoidal formula. 

An integral equation of the Volterra type is an equation which 
involves integrals of an unknown function x(t)/ which is a real or 
complex valueo function of a Single real variable t. The following 
IS a nonlinear integral equation of Volterra type; 

1 

Delves, L. M., "Numerical Software for Integral Equations," 

Nume rical Software Needs and Availability , Academic Press, pg. 305,119781. 
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x(t ) 


f(t) + 


1 . 

f* g ( t r S , X ( 


s) ) ds , where t > a. 


( 1 ) 


Functions f ano p are known, and the lower limit of intearation 

a oenctes a fixea real number while the unknown in this equation is 

the function x(t). Notice that the solution at any point t depends 

upon the values of x(s) on the interval a < s < t only. 

Another form of the Volterra equation is the nonlinear Volterra 

integral equation of convolution tyoe: 

t 


X (t ) 


f (t ) 


. f 


h(t-s/x(s))ds. 


( 2 ) 


0 


The linear form of this equation 

t 


X (t ) = f (t 3 


+ /'k (t , 


s 1 X ( s ) as 


(3) 


contains a function k called the kernel of the eouation. 

The equation, which is both linear ano of convolution type, 

t 


X Ct 1 


f (t ) 


, / 


a(t-slx(s)ds 




will be the focus of this paper. Function f(t) is called the forcing 
function. 

This linear Volterra equation of convolution type, which arises in 
a variety of mathematical, engineering and scientific applications, 

IS commonly referred to as the classic renewal equation. The forcing 
function f IS continuous, and the kernal a(tj is of class on each 
finite Subinterval of (0, OO). The renewal equation is defined as the 
expected number of renewals by time t. It can be usedf for example^ 
to oetermine the expected number of failures in a given time, or the 
probability distribution of the time at which all failures occur. 
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Due to tne strong link between Volterra equations and initial 
ue ordinary differential equation problems, the numerical methods 

equation (. 1 ) 

(5) 

relationshiD 

solution of 
the multi- 

step (predictor-corrector) methods. For an excellent comparison of 
existing algorithms see Delves in Jacobs C19781. 

SUFF VCLTEr^RA IMED^A L EQ U ATION .-- ITb V ULUI 1 0 f» KERNAL 

The type of Volterra equations for wnicn no algorithms exist are 
termed "stiff" Volterra problems. "Stiff equations contain time con- 
stants with greatly different values. Stability requires that the 
step size be not much greater than the smallest time constant present 
when conventional methods are usea. The terms due to the small time 

constants decay rapidly> but the step size cannot be increased because 

1 

ct staoility." 


usee to solve both are also strongly linxea. Therefore 
can be rewritten as 


X ( t + t ) 

0 


f (t + t ) 

0 




g(t+t »s,x(s))os 
0 


t tt 

r ® 

+ J g(t+t fSfX(s))os 

t 0 

which involves a step-by-step solution. This recursive 
contains implicitly the initial condition x(a) = f(a). 

Several numerical methods algorithms exist for the 
(1). The methods are often based on the Runge-Kutta or 


1 

Gear/ C.A.» "The Automatic Integration of Stiff Ordinary Differ- 
ential Equations''^ Informati on P rocessing 6o , pg. 187ftl9b8]. 
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example of the stiff Volterra equation 
exponential type of the form 

0 

Present at ion: 


t f 


max 


function graphically represented by 



larger than tf . Determine an 
max 

inction by choosing a truncation value 

ned over the interval [0/tf ]/ 

max 

the interval [Ofta ]. In order 

max 

olved in solving x(t3, given f(t) 
following graphs are presented. 


< c ompu t 1 ng x ( t ) 

n 


n 
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The following five points are illustrated by the previous graphs; 
1) Initially x(0) = f(0) = 1.0. 

?1 a(0) = Of therefore it is not necessary to know the value 

of x(t ) when solving for x(t ). 
n n 

3) For each successive computation of x(t where 0<t < tx 

n n ma X 

convoluted kernel a moves farther out on the time axis. 

ill tx is determined when either the function x equals tne 
max 

specified truncation value t or function x reaches a 

r 

steady state value. 

51 tx may be magnitudes larger than either tf or ta 

max max max 


In order to solve this "stiff" Volterra equation with convolution 

kernalf the step size must be carefully chosen to ensure stability. 

Therefore the smallest step size required by the functions involved 

in the integration determines the step size. But if this is used 

throuohout the computation of x(t)f the time required to compute the 

numoer of points necessary to define function x, from time zero to 

time tx f would not be computationally feasible due to the extremes 
max 

of the time maximums of each function involved in the computation. 
Storage of function xCt) could also cause a problem due to the re- 
cursive nature of the solution. 
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To illustrate the enormous number of calculations reauired to 


compute x(tl under these stiff con'd itions/ let tf = 1 minute and 

max 

ta = 1000 hours. Let the step size of function f equal 1 second, 
max 

yieldino 60 computed values for f(tj over the interval [0,tf ]. 

max 

If the 1 secona steo size were also used for function a, which would 

not be necessary, 3/o00,000 points would be computeo for function a. 

Instead, let the step size for function a equal 10 hours, yielding 

100 computed values for a(t) over the interval [0,ta ]. In oroer 

max 

to preserve stability, function x's step size must initially equal 1 
second. If this step size were used to compute x(tl over the inter- 
val 10, tx ], the number of computed x(t) values might exceed 3,600,000 
max 

values, due to the convoluted kernal. This methoa is not computation- 
ally feasiole because each new computeo x value requires the multi- 
plication of all Previous computed x values times the shifted function 
a. This woulo result in a minimum of 3,600,000 factorial multiplications 
plus 3,600,000 additions. 

One possible solution to the step size proolem woulo allow for 
the step size of function x(t) to increase as the rate of change 
of x(t) decreases by some specified amount. Doubling the step 
size each time an increase is warrantee would minimize the number of 
time values that must be stored with the computeo function values. 


IRAPLZ£Ii2iL_^LiiIi2N 


To obtain the approximate solution 
generalized trapezoidal formula, divide 
parts to yield the solution 

x(t I f(t ) + h/2 [x(0)aCt ) 


to equation (^), using the 
the interval [0,tJ into n 

+ 2(x[t )a(t -t ) +...+ 


n 


n n 1 n 1 

x(t )a(t )) + x(t )a(0n , 
n- 1 1 n 

where t = 0, t = h, ..., t = nh = t . 
0 1 n 
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x(t ) trust be eliminated from the right hand side of the equation 
n 

because this is the value currently being computea. Subtract the final 
proouct from both sides of the equation to yield 


x(t ) - h/2 [x(t )a(0}l f(t ) + h/d lxC0)a(t ) + 

n n n n 

2Cx(t )a(t -t ) +...+ 

1 n 1 

x(t )a(t ))] . 

n-1 1 

Factor out x(t ) to yield 


x(t )(1 - (h/2)a(0n ^ f(t ) + h/2 (x(O)aCt ) + 

n n n 

2(x(t )a(t “t ) +...+ 

1 n 1 

X (t )a(t ) )] . 

n-1 1 


Finally divide by (1 - (h/2)a(0)) to yield 


x(t - (h/2)a(0))] [f(t ) + h/2 (x(O)a(t ) + 

n n n 

2(x(t )a(t -t ) +...t 
1 n 1 

X Ct )a(t ))] . 
n-1 1 


This equation can be written in a more general form as 


x(t )::^[1/(1 - (h/2)a(0))J (f(t ) + (h/2)x(0)a(t ) + 

n n n 




X ( t ) a ( t -t ) ] 
1 n 1 


( 6 ) 


The following solution, based on the trapezoidal method described 
above» is a proposed solution to the stiff Volterra integral equation 
with convolution kernal: 
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11/(1 


X ( t ) 

n 


- ( (t -t )/21a(0) )] 

n n- 1 


IfCt ) + (t /2)x(0)a(t ) + 
n 1 n 


E 


( (t -t ) /a) X ( t ) • 
1+1 1-1 1 

a(t -t )] (7) 


n 1 

Notice that this solution ooes not require equal step sizes for 
functions X ana a; this is mandatory due to the extremes in the 
times for which the functions exist. The following is a graphical 
representation of the times involved in the previous solution for 
which functions a and x are known: 


I 


t 

k+1 
_l 


t 

n- 1 
_l 


a 


I I. 


time t 


t 

n 


X 


Function x has been computed up to time t and is currently 

k 

being solved for time t . In solution (7) a(t -t ) must be multi- 

n n n- 1 

plied by x(t ) in order to preserve stability and likewise 

n-1 

a(t -t 1 must be multiplied by x(t ) etc. It should be evident 
n n-2 n-2 

that a problem exists because x(t ) is the last computed value for 

k 

function where t < t < ... < t < t < t . Function x 

k k_tJ nj- 2 n-1 n 

coulo be extrapolated butr due to the recursive nature of the problem, 
this would introduce an inordinate amount of error, and instability 
woulo be introduced. 
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The following solution eliminates the ex t hapo 1 a t i on problem: 


x(t ) ifct ) + (t /2)x(0)a(t ) + 

n n 1 n 

k 

E ((t -t )/2)x(t )a(t -t ) + 

1=1 1+1 1-1 1 ni 

n- 1 

x(t ) • m ((t -t )/2)[(t -t 3/(t -t }]a(t -t n 

ki=k+l 1+1 1-1 ni nk ni 


n- 1 

1 - E (Ct -t l/2)[t -t )/t -t )]a(t -t ) - 
i=k+l 1+1 1-1 ik nk ni 

(Ct -t )/2)a(0) (8) 

n n- 1 


Solution (61 IS an expanded version of (. 7 ') } the integration has been 
split into two sections: 

1) The integration over the interval lO^t ], where function x 

k 

has been previously computed and can be interpolated if necessary; 

2 ) the integration over the interval It ,t 1/ where function x 

k n 

IS unknown, except for point x(t ). 

k 

This solution uses linear interpolation between the last known value 

x(t ) and unknown value x(t ) using the weights [(t - t )/(t -t }] and 
k n n 1 n k 

((t -t )/(t -t )1. Then x(t ) is eliminated from the right hand side of 
1 k n k n 

the equation in the exact same manner as described above for generating 
solution (bl. 

This solution minimizes the problems that make the stiff Volterra 
equation infeasiple to compute. It allows for the step size of function 
X to be initially as small as necessary to ensure stability, and then 
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enables the step size to increase as the function changes less oramati- 
cally. Each successive computed x(tj value uses all of the points of 
functions a and i, anu possibly other interpolated points of these 
functions, no matter how small the step sizes happen to be in relation 
to the current step size of function x. This method eliminates unnec- 
essary x(t) computations but incorporates all necessary points of 
functions a and f in each x(t) calculation. 

C ONCLUSIOf; 

No general purpose algorithms exist for solving the stiff Volterra 
integral eauation with convoluted kernal due partly to the problems of 
oetermininq the correct step size to ensure stability. Also, the con- 
voluted kernal requires that for each successive x(t) calculation 
all previous x values must be multiplied by the kernal function 

shifteo/ so that, in general^ a(t -t ) is multiplied by x(t ), 

n n-m n-m 

where t <t. 

n-m n 

The solution based on the trapezoidal method^ which allows for 
adaptive step sizes, eliminates these problems ana works well in 
solving the renewal equation with continuous, "well behaved" functions 
that are piece-wise 1 inear. 

Special thanks must be extended to Dr. J.J. Stiffler for his help 
with the final trapezoiaal solution. 
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